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Abstract 

Perturbations of stars and black holes have been one of the main topics of 
relativistic astrophysics for the last few decades. They are of particular im- 
portance today, because of their relevance to gravitational wave astronomy. In 
this review we present the theory of quasi-normal modes of compact objects 
from both the mathematical and astrophysical points of view. The discussion 
includes perturbations of black holes (Schwarzschild, Reissner-Nordstrom, Kerr 
and Kerr-Newman) and relativistic stars (non-rotating and slowly-rotating). 
The properties of the various families of quasi-normal modes are described, and 
numerical techniques for calculating quasi-normal modes reviewed. The suc- 
cesses, as well as the limits, of perturbation theory are presented, and its role 
in the emerging era of numerical relativity and supercomputers is discussed. 
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1 Introduction 



Helioseismology and asteroseismology are well known terms in classical astro- 
physics. From the beginning of the century the variability of Cepheids has been 
used for the accurate measurement of cosmic distances, while the variability of 
a number of stellar objects (RR Lyrae, Mira) has been associated with stel- 
lar oscillations. Observations of solar oscillations (with thousands of nonradial 
modes) hav e als o revealed a wealth of information about the internal structure 
of the Sun |204|. Practically every stellar object oscillates radially or nonradi- 
ally, and although there is great difficulty in observing such oscillations there 
are already results for various types of stars (O, B, ...). All these types of 
pulsations of normal main sequence stars can be studied via Newtonian theory 
and they are of no importance for the forthcoming era of gravitational wave 
astronomy. The gravitational waves emitted by these stars are extremely weak 
and have very low frequencies (cf. for a discussion of the sun [ [70| , and an im- 
portant new measurement of the sun's quadrupole moment and its application 
in the measurement of the anomalous precession of Mercury's perihelion |163|). 
This is not the case when wc consider very compact stellar objects i.e. neutron 
stars and black holes. Their oscillations, produced mainly during the formation 
phase, can be strong enough to be detected by the gravitational wave detectors 
(LIGO, VIRGO, GEO600, SPHERE) which are under construction. 

In the framework of general relativity (GR) quasi-normal modes (QNM) 
arise, as perturbations (electromagnetic or gravitational) of stellar or black hole 
spacetimes. Due to the emission of gravitational waves there are no normal 
mode oscillations but instead the frequencies become "quasi- normal" (complex) , 
with the real part representing the actual frequency of the oscillation and the 
imaginary part representing the damping. 

In this review we shall discuss the oscillations of neutron stars and black 
holes. The natural way to study these oscillations is by considering the linearized 
Einstein equations. Nev e rthe l ess, t here has been recent work on nonlinear black 
hole perturbations [101, 102, 103 104 , 10C| while, as yet nothing is known for 
nonlinear stellar oscillations in general relativity. 

The study of black hole perturbations was initiated by the pioneering work 
of Regge and Wheeler |173] in the late 50s and was continued by Zerilli |212|. 
The perturbations of relativistic st ars in GR wer e firs t studied in the late 60s by 
Kip Thorne and his coUaborators (20| |l9|, |l9|, |20§. The initial aim of Regge 



and Wheeler was to study the stability of a black hole to small perturbations 
and they did not try to connect these perturbations to astrophysics. In con- 
trast, for the case of relativistic stars, Thome's aim was to extend the known 
properties of Newtonian oscillation theory to general relativity, and to estimate 
the frequencies and the energy radiated as gravitational waves. 

QNMs were first pointed out by Vishveshwara |207| in calculations of the 
scattering of gravitational waves by a Schwarzschild black hole, while Press [164| 
coined the term quasi-normal frequencies. QNM oscillations have been found 
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in perturbation calculations of particles falling into Schwarzschild [ J73| a nd Kerr 
black holes Q and in the collapse of a star to form a black hole |]66|, |67| , |68t . 
Numerical investigations of the fully nonlinear equations of general relativity 
have provided results which agree with the results of perturbation calculations; 
in particular numerical studies of the head-on collision of t wo b lack holes 1 30 , 2^ 
(cf. Figure |l|) and gravitational c ollap se to a Kerr hole |191|. Recently, Price, 
PuUin and collaborators |170, 101, ^ have pushed forward the agreement 



between full nonlinear numerical results and results from perturbation theory 
for the collision of two black holes. This proves the power of the perturbation 
approach even in highly nonlinear problems while at the same time indicating 
its limits. 

In the concluding remarks of their pioneering paper on nonradial oscillations 
of neutron stars Thorne and CampoUataro [202| described it as "jwsf a modest 
introduction to a story which promises to he long, complicated and fascinating^^ . 
The story has undoubtedly proved to be intriguing, and many authors have 
contributed to our present understanding of the pulsations of both black holes 
and neutron stars. Thirty years after these prophetic words by Thorne and 
CampoUataro hundreds of papers have been written in an attempt to understand 
the stability, the characteristic frequencies and the mechanisms of excitation of 
these oscillations. Their relevance to the emission of gravitational waves was 
always the basic underlying reason of each study. An account of all this work 
will be attempted in the next sections hoping that the interested reader will 
find this review useful both as a guide to the literature and as an inspiration for 
future work on the open problems of the field. In the next section we attempt to 
give a mathematical definition of QNMs. The third and fourth section will be 
devoted to the study of the black hole and stellar QNMs. In the fifth section we 
discuss the excitation and observation of QNMs and finally in the sixth section 
we will mention the more significant numerical techniques used in the study of 
QNMs. 
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Figure 1: QNM ringing after the head-on collision of two unequal mas s black 
holes l2^. The continuous line corresponds to the full nonlinear numerical 
calculation while the dotted line is a fit to the fundamental and first overtone 
QNM. 

2 Normal Modes — Quasi-Normal Modes — Res- 
onances 

Before discussing quasi-normal modes it is useful to remember what normal 
modes are! 

Compact classical linear oscillating systems such as finite strings, mem- 
branes, or cavities filled with electromagnetic radiation have preferred time 
harmonic states of motion (lu is real): 



if dissipation is neglected. (We assume x to be some complex valued field.) 
There is generally an infinite collection of such periodic solutions, and the "gen- 
eral solution" can be expressed as a superposition, 



of such normal modes. The simplest example is a string of length L which 
is fixed at its ends. All such systems can be described by systems of partial 
differential equations of the type (x may be a vector) 



Xn{t,x) =e"^"'Xn{x), 1,2,3... 



(1) 



OO 




(2) 



n=l 



dt 



(3) 
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where A is a linear operator acting only on the spatial variables. Because 
of the finiteness of the system the time evolution is only determined if some 
boundary conditions are prescribed. The search for solutions periodic in time 
leads to a boundary value problem in the spatial variables. In simple cases it 
is of the Sturm-Liouville type. The treatment of such boundary value problems 
for differential equations played an important role in the development of Hilbert 
space techniques. 

A Hilbert space is chosen such that the differential operator becomes sym- 
metric. Due to the boundary conditions dictated by the physical problem, A 
becomes a self-adjoint operator on the appropriate Hilbert space and has a pure 
point spectrum. The eigenfunctions and eigenvalues determine the periodic 
solutions (|l|). 

The definition of self-adjointness is rather subtle from a physicist's point of 
view since fairly complicated "domain issues" play an essential role. (See 
where a mathematical exposition for physicists is given.) The wave equation 
modeling the finite string has solutions of various degrees of differentiability. 
To describe all "realistic situations", clearly C°° functions should be sufficient. 
Sometimes it may, however, also be convenient to consider more general solu- 
tions. 

From the mathematical point of view the collection of all smooth functions is 
not a natural setting to study the wave equation because sequences of solutions 
exist which converge to non-smooth solutions. To establish such powerful state- 
ments like (||) one has to study the equation on certain subsets of the Hilbert 
space of square integrable functions. For "nice" equations it usually happens 
that the eigenfunctions are in fact analytic. They can then be used to gen- 
erate, for example, all smooth solutions by a pointwise converging series (^. 
The key point is that we need some mathematical sophistication to obtain the 
"completeness property" of the eigenfunctions. 

This picture of "normal modes" changes when we consider "open systems" 
which can lose energy to infinity. The simplest case are waves on an infinite 
string. The general solution of this problem is 

X{t,x)^ A{t-x) + B{t + x) (4) 

with "arbitrary" functions A and B. Which solutions should we study? Since 
we have all solutions, this is not a serious question. In more general cases, 
however, in which the general solution is not known, wc have to select a certain 
class of solutions which we consider as relevant for the physical problem. 

Let us consider for the following discussion, as an example, a wave equation 
with a potential on the real line, 

^X+(-^ + nx))x^O. (5) 

Cauchy data x(0, x), 9tx(0, x) which have two derivatives determine a unique 
twice differentiable solution. No boundary condition is needed at infinity to 
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determine the time evolution of the data! This can be estabHshed by fairly 



simple PDE theory |116 



There exist solutions for which the support of the fields are spatially compact, 
or - the other extreme - solutions with infinite total energy for which the fields 
grow at spatial infinity in a quite arbitrary way! 

From the point of view of physics smooth solutions with spatially compact 
support should be the relevant class - who cares what happens near infinity! 
Again it turns out that mathematically it is more convenient to study all solu- 
tions of finite total energy. Then the relevant operator is again self-adjoint, but 
now its spectrum is purely "continuous" . There are no eigenfunctions which are 
square integrable. Only "improper eigenfunctions" like plane waves exist. This 
expresses the fact that we find a solution of the form (|^) for any real uj and 
by forming appropriate superpositions one can construct solutions which are 
"almost eigenfunctions". (In the case V{x) = these are wave packets formed 
from plane waves.) These solutions are the analogs of normal modes for infinite 
systems. 

Let us now turn to the discussion of "quasi-normal modes" which are concep- 
tually different to normal modes. To define quasi-normal modes let us consider 
the wave equation for potentials with V >Q which vanish for \x\ > xq. Then 
in this case all solutions determined by data of compact support are bounded: 
|x(i, < C We can use Laplace transformation techniques to represent such 
solutions. The Laplace transform x{s,x) (s > real) of a solution a:) is 

POO 

X{s,x) - / e-''x{t,x)dt, (6) 
Jo 

and satisfies the ordinary differential equation 



"'x - x" + Vx^ +sx(0, x) + dtxio, x), (7) 



s 



where 



s'x -x' + Vx^Q (8) 



is the homogeneous equation. The boundedncss of x implies that x is analytic 
for positive, real s, and has an analytic continuation onto the complex half plane 
Re{s) > 0. 

Which solution x of this inhomogeneous equation gives the unique solution 
in spacetime determined by the data? There is no arbitrariness; only one of the 
Green functions for the inhomogeneous equation is correct! 

All Green functions can be constructed by the following well known method. 
Choose any two linearly independent solutions of the homogeneous equation 
f-{s,x) and /+(s,a;), and define 

^ - ^ f-{s,x')f+{s,x) (x' < x), . , 

T4^(s) A(s,a;)/+(s,a;') {x' > x), 
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where W{s) is the Wronskian of /_ and /+. If we denote the inhomogeneity 
of by j, a solution of (0) is 

/•oo 

X{s,x)^ G{s,x,x)j{s,x')dx' . (10) 



We stih have to select a unique pair of solutions f-,f+. Here the information 
that the solution in spacetime is bounded can be used. The definition of the 
Laplace transform implies that x is bounded as a function of x. Because the 
potential V vanishes for > xq, the solutions of the homogeneous equation (||) 
for |a;| > xq are 

/ = e±^^ (11) 

The following pair of solutions 

/+ = e"'*^ ior x > Xq, /_ = 6+"*^ for x < -a;o, (12) 

which is linearly independent for Re{s) > 0, gives the unique Green fmiction 
which defines a bounded solution for j of compact support. Note that for 
Re{s) > the solution /+ is exponentially decaying for large x and /_ is expo- 
nentially decaying for small x. For small x however, /+ will be a linear com- 
bination a(s)e~*^ -1- b{s)e^'^ which will in general grow exponentially. Similar 
behavior is found for 

Quasi-Normal mode frequencies s„ can be defined as those complex numbers 
for which 

/+(s„,a;) = c(s„)/_(s„,a;), (13) 

that is the two functions become linearly dependent, the Wronskian vanishes 
and the Green function is singular! The corresponding solutions /+(s„,x) are 
called quasi eigenfunctions. 

Are there such numbers s„? From the boundedness of the solution in space- 
time we know that the unique Green function must exist for Re{s) > 0. Hence 
/+,/_ are linearly independent for those values of s. However, as solutions 
/+ , /_ of the homogeneous equation (|^) they have a unique continuation to the 
complex s plane. In p5[ it is shown that for positive potentials with compact 
support there is always a countable number of zeros of the Wronskian with 
Re{s) < 0. 

What is the mathematical and physical significance of the quasi-normal fre- 
quencies Sn and the corresponding quasi- normal functions /+? First of all we 
should note that because of Re{s) < the function /+ grows exponentially 
for small and large x\ The corresponding spacetime solution e*"*/+(sn, a;) is 
therefore not a physically relevant solution, unlike the normal modes. 

If one studies the inverse Laplace transformation and expresses x as a com- 
plex hne integral (a > 0), 

1 /■+°° 

X{t, x)^— e^^+^'^'xia + IS, x)ds, (14) 

27ri ,1^, 
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one can deform the path of the complex integration and show that the late time 
behavior of solutions can be approximated in finite parts of the space by a finite 
sum of the form 



Here we assume that i?e(s„+i) < i?e(s„) < 0, s„ = a„ + i/3„. The approxi- 
mation ~ means that if we choose Xq, Xi, e and to then there exists a constant 
C(to,Xo,Xi,e) such that 



holds for t > to, Xq < X < xi, e > with C{to,xo,xi,e) independent of t. 
The constants a„ depend only on the data jSSt ! This implies in particular that 
all solutions defined by data of compact support decay exponentially in time 
on spatially bounded regions. The generic leading order decay is determined 
by the quasi- normal mode frequency with the largest real part si, i.e. slowest 
damping. On finite intervals and for late times the solution is approximated by 
a finite sum of quasi eigenfunctions (|l5|) . 

It is presently unclear whether one can strengthen ( ]l6| ) to a statement 
like (|^), a pointwise expansion of the late time solution in terms of quasi- normal 
modes. For one particular potential (Poschl- Teller) this has been shown by 
Beyer [||. 

Let us now consider the case where the potential is positive for all x, but 
decays near infinity as happens for example for the wave equation on the static 
Schwarzschild spacet ime. Data of compact support determine again solutions 
which are bounded |ll7|l . Hence we can proceed as before. The first new 
point concerns the definitions of f±. It can be shown that the homogeneous 
equation (^) has for each real positive s a unique solution f^(s,x) such that 
lima;^oo(e*^/+(s, x)) = 1 holds and correspondingly for These functions 
are uniquely determined, define the correct Green function and have analytic 
continuations onto the complex half plane Re{s) > 0. 

It is however quite complicated to get a good representation of these func- 
tions. If the point at infinity is not a regular singular point, we do not even get 
converging series expansions for f± . (This is particularly serious for values of s 
with negative real part because we expect exponential growth in x). 

The next new feature is that the analyticity properties of f± in the complex 
s plane depend on the decay of the potential. To obtain information about 
analytic continuation, even use of analyticity properties of the potential in x is 
made! Branch cuts may occur. Nevertheless in a lot of cases an infinite number 
of quasi-normal mode frequencies exists. 

The fact that the potential never vanishes may, however, destroy the expo- 
nential decay in time of the solutions and therefore the essential properties of 



N 




(15) 



n=l 




(16) 
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the quasi- normal modes. This probably happens if the potential decays slower 
than exponentially. There is, however, the following way out: Suppose you want 
to study a solution determined by data of compact support from i = to some 
large finite time t = T. Up to this time the solution is - because of domain of 
dependence properties - completely independent of the potential for sufficiently 
large x. Hence we may see an exponential decay of the form ( p^ in a time 
range ti < t < T . This is the behavior seen in numerical calculations. The 
situation is similar in the case of a-decay in quantum mechanics. A comparison 
of quasi-normal modes of wave equations and resonances in quantum theory can 
be found in the appendix, see section |^. 
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3 Quasi-Normal Modes of Black Holes 



One of the most interesting aspects of gravitational wave detection will be the 
connection with the existence of black holes |201]. Although there are presently 
several indirect ways of identifying black holes in the universe, gravitational 
waves emitted by an oscillating black hole will carry a unique fingerprint which 
would lead to the direct identification of their existence. 

As we mentioned earlier, gravitational radiation from black hole oscillations 
exhibits certain characteristic frequencies which are independent of the pro- 
cesses giving rise to these oscillations. These "quasi-normal" frequencies are 
directly connected to the parameters of the black hole (mass, charge and angu- 
lar momentum) and for stellar mass black holes are expected to be inside the 
bandwidth of the constructed gravitational wave detectors. 

The perturbations of a Schwarzschild black hole reduce to a simple wave 
equation which has been studied extensively. The wave equation for the case 
of a Rcissner-Nordstrom black hole is more or less similar to the Schwarzschild 
case, but for Kerr one has to solve a system of coupled wave equations (one for 
the radial part and one for the angular part). For this reason the Kerr case has 
been studied less thoroughly. Finally, in the case of Kerr-Newman black holes 
we face the problem that the perturbations cannot be separated in their angular 
and radial parts and thus apart from special cases \124\ the problem has not 
been studied at all. 



3.1 Schwarzschild Black Holes 

The study of perturbations of Schwarzschild black holes assumes a small per- 
turbation hf^i, on a static spherically symmetric background metric 

ds^ = g^^^dx^'dx" = -e^'^^'^dt^ + e^^^'^dr^ + [dO^ + sin^ ed(j)'^) , (17) 

with the perturbed metric having the form 

9iiu ^ + (18) 

which leads to a variation of the Einstein equations i.e. 

6G^,,^Att6T^,. (19) 

By assuming a decomposition into tensor spherical harmonics for each h^^ of 
the form 

X(t, r, e,^)^Y. ^^^^^Y^mie. 0), (20) 

the perturbation problem is reduced to a single wave equation, for the func- 
tion XimifTt) (which is a combination of the various components of h^^). It 
should be pointed out that equation (pQ) is an expansion for scalar quantities 
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only. From the 10 independent components of the h^^i, only hu, htr, and hrr 
transform as scalars under rotations. The hte, ht^, hre, and hrrf, transform as 
components of two-vectors under rotations and can be expanded in a series of 
vector spherical harmonics while the components hgg, hstf,, and h^^ transform 
as components of a 2x2 tensor an d can be expanded in a series of tensor spher- 
ical harmonics (see |202| , |212[ |152[ for details). There are two classes of vector 
spherical harmonics (polar and ELxial) which are build out of combinations of 
the Levi-Civita volume form and the gradient operator acting on the scalar 
spherical harmonics. The difference between the two families is their parity. 
Under the parity operator tt a spherical harmonic with index £ transforms as 
(— 1)^, the polar class of perturbations transform under parity in the same way, 
as (—1)^, and the axial perturbations as (— 1)^^^0. Finally, since we are dealing 
with spherically symmetric spacetimes the solution will be independent of m, 
thus this subscript can be omitted. 

The radial component of a perturbation outside the event horizon satisfies 
the following wave equation, 

^X.+ (-^ + 14(r)jx. = 0, (21) 

where is the "tortoise" radial coordinate defined by 

r,^r + 2A/log(r/2Af - 1), (22) 

and M is the mass of the black hole. 
For "axial" perturbations 

.,M^(i-M) 

is the effe ctive potential or (as it is known in the literature) Regge- Wheeler 
potential [[l73| , which is a single potential barrier with a peak around r = 3M, 
which is the location of the unstable photon orbit. The form ( p3| ) is true even if 
we consider scalar or electromagnetic test fields as perturbations. The parameter 
a takes the values 1 for scalar perturbations, for electromagnetic perturbations, 

^In the literature the polar perturbations are also called even-parity because they are 
characterized by their behavior under parity operations as discussed earlier, and in the same 
way the axial perturbations are called odd-parity. We will stick to the polar/axial terminology 
since there is a confusion with the definition of the parity operation, the reason is that to 
most people, the words "even" and "odd" imply that a mode transforms under tt as (—1)-^" 
or (— respectively (for n some integer). However only the polar modes with even ( 
have even parity and only axial modes with even £ have odd parity. If £ is odd, then polar 
modes have odd parity and axial modes have even parity. Another terminology is to call the 
polar perturbations spheroidal and the axial ones toroidal. This definition is coming from the 
study of stellar pulsations in Newtonian theory and represents the type of fluid motions that 
each type of perturbation induces. Since we are dealing both with stars and black holes we 
will stick to the polar /axial terminology. 



£{£ + 1) 2aM 



(23) 
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and —3 for gravitational perturbations and can be expressed as cr — 1 — s^, where 
s = 0, 1, 2 is the spin of the perturbing field. 



For "polar" perturbations the effective potential was derived by Zerilli [212| 
and has the form 

''^^^=\^-—) .3(„, + 3M)2 ^ (24) 

where 

2n = l)(^ + 2). (25) 

Chandrasekhar has shown that one can transform the equation (|2^) for 
"axial" modes to the corresponding one for "polar" modes via a transforma- 
tion involving differential operations. It can also be shown that both forms 
are connected to the Bardeen-Press |Q perturbation equation derived via the 
Newman-Penrose formalism. The potential 14 (r*) decays exponentially near 
the horizon, r* — oo, and as for +00. 

From the form of equation ( pT|) it is evident that the study of black hole 
perturbations will follow the footsteps of the theory outlined in section |2[ 



Kay and Wald [117| have shown that solutions with data of compact sup- 
port are bounded. Hence we know that the time independent Green function 
G(^s^r^,^r'^') is analytic for Re{s) > 0. The essential difficulty is now to obtain 
the solutions f± (cf. equation (^0|)) of the equation 

s'x - X" + Vx = 0, (26) 

(prime denotes differentiation with respect to ) which satisfy for real, positive 
s: 

f+ - e-'''' for n ^00, /_ - e+''*"^ for r, -00. (27) 

To determine the quasi-normal modes we need the analytic continuations of 
these functions. 

As the horizon (r* — > 00) is a regular singular point of (|2^), a representation 
of /-(r*, s) as a converging series exists. For ill = i it reads: 

00 

/-(r,,s) = (r-ir^a„(s)(r-l)". (28) 

Tl = 

The series converges for all complex s and |r — 1| < 1 |162|. (The analytic 



extension of /_ is investigated in [ 1 1 5 1 . ) The result is that /_ has an extension to 
the complex s plane with poles only at negative real integers. The representation 
of /+ is more complicated: Because infinity is a singular point no power series 
expansion like (Eq) exists. A representation coming from the iteration of the 



defining integral equation is given by Jensen and Candelas [ |115| , see also [159|. 
It turns out that the continuation of /+ has a branch cut Re{s) < due to the 

„-2 



decay r ^ for large r [115| 
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The most extensive mathematical investigation of quasi-normal modes of 
the Schwarzschild solution is contained in the paper by Bachclot and Motet- 
Bachelot [^^. Here the existence of an infinite number of quasi- normal modes 
is demonstrated. Truncating the potential ( p^ ) to make it of compact support 
leads to the estimate (p^. 

The decay of solutions in time is not exponential because of the weak decay 
of the potential for large r. At late ti mes, the quasi-normal oscillations are 
swamped by the radiative tail [166, 167 1. This tail radiation is of interest in its 
own right since it originates on the background spacetime. The first authorita- 
tive study of nearly spherical collapse, exhibiting radiative tails, was performed 
by Price |l6|, [l67| . 

Studying the behavior of a massless scalar field propagating on a fixed 
Schwarzschild background, he showed that the field dies off with the power- 
law tail, 

(29) 

at late times, where P = 1 if the field is initially static, and P — 2 otherwise. 
This behavior has been seen in various calculations, for example the gravita- 
tional collapse simulations by Cunningham, Price and Moncrief |6^, ^ 
Today it is apparent in any simulation involving evolutions of various fields 
on a black hole background including Schwarzschild, Reissner- Nordstrom [106|, 
and Kerr [132, 133|. It has also been observed in simulations of axial oscil- 
lations of ne utron stars [|lj], and should also be present for polar oscillations. 
Leaver [ 136 has studied in detail these tails and associated this power low tail 
with the branch-cut integral along the negative imaginary u axis in the com- 
plex CO plane. His suggestion that there will be radiative ta ils o bservable at J'^ 
and 7i+ has been verified by Gundlach, Price, and PuUin [106|. Similar results 
were arrived at recently by Ching et al. [|62[ in a more extensive study of the 
late time behavior. In a nonlinear study Gundlach, Price, and PuUin [107[ have 
shown that tails develop even when the collapsing field fails to produce a black 
hole. Finally, for a study of tails in the presence of a cosmological constant refer 
to [^, while for a recent study, using analytic methods, of the late-time tails of 
linear scalar fields outside Schwarzschild and Kerr black holes refer to |36, ^7\ . 



Using the properties of the waves at the horizon and infinity given in equa- 
tion ( p7| ) one can search for the quasi-normal mode frequencies since practically 
the whole problem has been reduced to a boundary value problem with s = iu) 
being the complex eigenvalue. The procedure and techniques used to solve the 
problem will be discussed later in section |^, but it is worth mentioning here 
a simple approach to calculate the QNM frequencies proposed by Schutz and 
Will [18C[. The approach is based on the standard WKB treatment of wave 
scattering on the peak of the potential barrier, and it can be easily shown that 
the complex frequency can be estimated from the relation 

1/2 



(" + 









(30) 
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where ro is the peak of the potential barrier. For £ = 2 and n = (the funda- 
mental mode) the complex frequency is Mu) k, (0.37, —0.09), which for a IOM0 
black hole corresponds to a frequency of 1.2 kHz and damping time of 0.55 ms. 
A few more QNM frequencies for £ = 2,3 and 4 are listed in table |l|. 



n 


£ = 2 


£ = i 


£ = A 



1 
2 
3 


0.37367 -0.08896 i 
0.34671 -0.27391 i 
0.30105 -0.47828 i 
0.25150 -0.70514 i 


0.59944 -0.09270 i 
0.58264 -0.28130 i 
0.55168 -0.47909 i 
0.51196 -0.69034 i 


0.80918 -0.09416 i 
0.79663 -0.28443 i 
0.77271 -0.47991 i 
0.73984 -0.68392 i 



Table 1: The first four QNM frequencies (toM) of the Schwarzschild black hole 
for £ = 2,3, and 4 Jj,?^ /- The frequencies are given in geometrical units and for 
conversion into kHz one should multiply by 27r(5142iJz) x {Mq/M). 

Figure H shows some of the modes of the Schwarzschild black hole. The num- 
ber of modes for each harmonic index t is infinite, as was mathematically proven 
by Bachelot and Motet-Bachelot |^^. This was also implied in an earlier work 
by Ferrari and Mashhoon psf , and it has been seen in the numerical calculations 
It can be also seen that the imaginary part of the frequency grows 



157| 



very quickly. This means that the higher modes do not contribute significantly 
in the emitted gravitational wave signal, and this is also true for the higher £ 
modes (octapole etc.). As is apparent in figure || that there is a special purely 
imaginary QNM frequency. The existence of "algebraically special" solutions 
for perturbations of Schwarzschild, Reissner-Nordstrom and Kerr black holes 
were first pointed out by Chandrasekhar jSTj. I t is s till questionable whether 
these frequencies should be considered as QNMs 137 and there is a suggestion 
that the potential might become tran sparent for these frequencies |l^ . For a 
more detailed discussion refer to |144|. 

As a final comment we should mention that as the order of the modes in- 
creases the real part of the frequency remains constant, while t he im aginary 
part increases proportionally to the order of the mode. NoUert [157| derived 
the following approximate formula for the asymptotic behavior of QNMs of a 
Schwarzschild black hole, 



MoJr, ~ 0.0437- 



71 



(2n + 1)1/2 



--{2n 



1) 



71 



(2n + 1)1/2 



, (31) 



where 71 = 0.343, 0.7545 and 2.81 for £ = 2, 3 and 6, corr espondingly, and 
n ^ 00. The above relation was later verified in jl^ and [ 143 1 . 



For large values of £ the distribution of QNMs is given by [164, g6| g5| 113 

' £-\ i[n + - 

2 \ 2 



3V3Mu 



(32) 
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Figure 2: The spectrum of QNM for a Schwarzschild black-hole, for 1 = 2 
(diamonds) and £ = 3 (crosses) /p^. The 9th mode for t = 2 and the 41st for 
i = 3 are "special", i.e. the real part of the frequency is zero (s = ito). 



For a mathematical proof refer to . 

The perturbations of Reissner-Nordstrom black holes, due to the spherical 
symmetry of the solution, follow the footsteps of the analysis that we have 
presented in this section. Most of the work was done during the seventies by 



Zerilli [213|, Moncrief |153, 154| and later by Chandrasekhar and Xanthopou- 
los 20£]. For an extensive discussion refer to |Q. We have again wave equa- 
tions of the form (21), one for each parity with potentials which are like (P3| 



and (24) plus extra terms which relate to the charge of the black hole. An inter- 
esting feature of the charged black holes is that any perturbation of the gravi- 
tational (electromagnetic) field will also induce electromagnetic (gravitational) 
perturbations. In other words, any perturbation of the Reissner-Nordstrom 
spacetime will produce both electromagnetic and gravitational radiation. Again 
it has been shown that the solutions for the odd parity oscillations can be de- 
duced from the solutions for even parity oscillations and vice versa |55 . The 



QNM frequencies of the Reissner-Nordstrom black hole have been calculated 
by Gunter ]10^, Kokkotas, and Schutz [1129|, Leaver |137ll, Andersson 1^, and 



Kokkotas, and Schutz [129|, Leaver |137|, Andersson 
lately for the nearly extreme case by Andersson and Onozawa p6||. 
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3.2 Kerr Black Holes 



The Kerr metric represents an axisymmetric, black hole solution to the source 
free Einstein equations. The metric in {t,r,6,ip) coordinates is 



1 



2Mr 



2Ma\ sin^ 6 



AMar sin^ 9 



-dtdip 



-dr^ + Y.de^ 



sin^ 6dip'^ 



with 



2Mr 



a? cos^ 9. 



(33) 



(34) 



M is the mass and Q < a < M the rotational parameter of the Kerr metric. 
The zeros of A are 

M±{M^ -a'Y'^, (35) 

and determine the horizons. For r_|_ < r < oo the spacetime admits locally a 
timelike KiUing vector. In the ergosphere region 



r+ <r < M 



cos^ I 



xl/2 



(36) 



the Killing vector d/dt which is timelike at infinity, becomes spacelike. The 
scalar wave equation for the Kerr metric is 



(r2+a2)2 



dr 



- sin^ ( 



A-+i^)- 
dr 



1 



AMar d'^x 



A dtdif 



d 



sin 9 89 



sin.|)-2. 



1_ 

a" ^ 

a{r~ M) 



d^X 



-2cr 



— r — la cos ( 



dx 

at 



fcr^ cot^ ( 



A 

-^)x = o, 



I cost 



sm 



dx 

dip 

(37) 



where a = 0,±1,±2 for scalar, electromagnetic or gravitational perturbations, 
respectively. As the Kerr metric outside the horizon (r > r^) is globally hy- 
perbolic, the Cauchy problem for the scalar wave equation ( |37| ) is well posed 
for data on any Cauchy surface. However, the coefficient of d^x/9'p'^ becomes 
negative in the ergosphere. This implies that the time independent equation we 
obtain after the Fourier or Laplace transformation is not elliptic! 

For linear hyperbolic equations with time independent coefficients, we know 
that solutions determined by data with compact support are bounded by ce'''*, 
where 7 is independent of the data. It is not known whether all such solutions 
are bounded in time, i.e. whether they are stable. 

Assuming harmonic time behavior x = e'"*x('': V'); ^ separation in angular 



(38) 



and radial variables was found by Teukolsky ]196| |: 

X{r,9,(p) = R(r,uj)S{9,uj)e'""^. 
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Note that in contrast to the case of spherical harmonics, the separation is 



dependent. To be a solution of the wave equation (37), the functions R and S 
must satisfy 



1 



sin e de 



dr 

where K = 



dO 



a^iJ^ cos^ Q + 2aa;cr cos Q 



m'^+a^+2ma cos 6 



+E 



^dR 
dr 



5 = 0, 
(39) 

2i(j{r - 1)K - A{4:iaru; + \)] R = 0, (40) 



am, 



X = E + a'^uj'^ + 2mnuj — a{a + 1), and E is the 



separation constant. 

For each complex a^cj^ and positive integer m, equation (39) together with 
the boundary conditions of regularity at the axis, determines a singular Sturm- 
Liouville eigenvalue problem. It has solutions for eigenvalues E{£,m^,a^uj^), 
\m\ < I. The eigenfunctions are the spheroidal (oblate) harmonics S'£|„j|(6'). 
They exist for all complex w^. For real uj^ the spheroidal harmonics are com- 
plete in the sense that any function of z = cos 0, absolutely intcgrable over the 



interval [—1,1], can be expanded into spheroidal harmonics of fixed m [181|. 
Furthermore, functions A{d^ Lp) absolutely integrable over the sphere can be 
expanded into 



+1 



irmp 



(41) 



1=0 m=- 



For general complex such an expansion is not possible. There is a co untable 
number of "exceptional values" where no such expansion exists |14^ . 

Let us pick one such solution S^^^ {0, uj) and consider some solutions R{r, uj) 
of dl^) with the corresponding E{e, |m|, w). Then (r, w)5f|„|e™'^e-''^* is a 
solution of (^7|). Is it possible to obtain "all" solutions by summing over £, m and 
integrating over w? For a solution in spacetime for which a Fourier transform 
in time exists at any space point (square integrable in time), we can expand 
the Fourier transform in spheroidal harmonics because lo is real. The coefficient 
i?(r, UJ, m) will solve equation (40). Unfortunately, we only know that a solution 
determined by data of compact support is exponentially bounded. Hence we 
can only perform a Laplace transformation. 

We proceed therefore as in section |2[ Let x{s,r,9,(p) be the Laplace trans- 
form of a solution determined by data x(i, J', ¥') = and dtx{t,r,0,ip) = p, 
while X is analytic in s for real s > 7 > 0, and has an analytic continuation 
onto the half-plane Re{s) > 7. For real s we can expand x i'^to a converging 
sum of spheroidal harmonics |14S] 



X{s,r,d,(p) 



00 +e 

EE 

e=o m=-t 



R{e,m,s)Sa,n\{ 



(42) 
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R{(, TO, s) satisfies the radial equation with iuj — s. This representation of 
X does, however, not hold for all complex values in the half-plane on which it is 
defined. Nevertheless is it true that for all values of s which arc not exceptional 
an expansion of the form exists^ 

To define quasi-normal modes we first have to define the correct Green 
function of ( ^o| ) which determines R{£, to, s) from the data p. As usual, this 
is done by prescribing decay for real s for two linearly independent solutions 
^R{i, TO, s) and analytic continuation. Out of the Green function for R{i, to, s) 
and S'^iml (— s^j ^) for real s we can build the Green function G{s, r, r' , 0, 9' , (j>, (j)') 
by a series representation like ( ^2[ ) . Analytic continuation defines G on the half 
plane Re{s) > 7. For non exceptional s we have a series representation. (We 
must define G by this complicated procedure because the partial differential 
operator corresponding to (^), ( ^0|) is not elliptic in the ergosphere.) 

Normal and quasi-normal modes appear as poles of the analytic continuation 
of G. Normal modes are determined by poles with Re{s) > and quasi- normal 
modes by Re{s) < 0. Suppose all such values are different from the exceptional 
values. Then we have always the series expansion of the Green function near 
the poles and we see that they appear as poles of the radial Green function of 
R{£, TO, s). 

To relate the modes to the asymptotic behavior in time we study the inverse 
Laplace transform and deform the integration path to include the contributions 
of the poles. The decay in time is dominated cither by the normal mode with the 
largest (real) eigenfrequency or the quasi-normal mode with the largest negative 
real part. 

It is apparent that the calculation of the QNM frequencies of the Kerr black 
hole is more involved than the Schwarzschild and Reissner-Nordstrom cases. 
This is the reason that there have only been a few attempts [135, 185, 123, 160| 
in this direction. 

The quasi-normal mode frequencies of the Kerr-Newman black hole have 
not yet been calculated, although they are more general than all other types 
of perturbation. The reason is the complexity of the perturbation equations 
and, in particular, their non-separability. This can be understood through the 
following analysis of the perturbation procedure. The equations governing a 
perturbing massless field of spin a can be written as a set of 2(7 -f 1 wavelike 
equations in which the various different helicity components of the perturbing 
field are coupled not only with each other but also with the curvature of the 
background space, all with four independent variables as coordinates over the 
manifold. The standard problem is to decouple the 2a + 1 equations or at least 
some physically important subset of them and then to separate the decoupled 
equations so as to obtain ordinary differential equations which can be handled 

■^The definition of S£|„|(— on the complex plane is made unique by fixing certain 
conventions about the branch cuts. The exceptional points are the beginnings of branch cuts. 
R{1, m, s) as a solution of ( ^o[ ) is defined also at the exceptional points; just the series does 
not exist there. 
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by one of the previously stated methods. For a discussion and estimation of the 
QNM frequencies in a restrictive case refer to [124|. 



3.3 Stability and Completeness of Quasi-Normal Modes 

From the normal modes one can learn a lot about stability. Take as an example 
linear stellar oscillation within the framework of Newton's theory of gravity. As 
outlined at the beginning of section ^ we have a sequence of normal modes 
ujn with real. The general solution is a convergent linear combination of 
the corresponding eigenfunctions. Hence all eigenfunctions are bounded if lu is 
purely imaginary, i.e. uj^ < 0. 

Therefore the spectrum contains all the information about stability. To 
discuss stability for systems with quasi-normal modes, let us consider a case 
like equation (j^) with the assumption that V is of compact support but not 
necessarily positive. 

Data of compact support define solutions which grow at most exponentially 
in time 

m,x)\<ce''\ (43) 

where a is independent of the data. As outlined in appendix ^ eigenvalues 
necessarily have s„ > and the eigenfunctions determine solutions growing 
exponentially in time. If no eigenvalues exist, the solution can not grow expo- 
nentially. Polynomial growth is still possible and related to the properties of the 
Laplace transform of the Green function at s = 0. As the potential has compact 
support, the functions /±(s, x) are analytic for all s. Hence, the Green fmiction 
can at most have a pole at s = 0. A pole of order two and higher implies poly- 
nomial growth in time. If the potential is positive, energy conservation shows 
that the field can grow at most linearly in time and therefore we can have at 
most a pole of order 2 at s = 0. 

If we define stability as boundedness in time for all solutions with data of 
compact support, properties of quasi-normal modes can not decide the stability 
issue. However, the appearance of a normal mode proves instability. If the 
support of the potential is not compact everything becomes more complicated. 
In particular, it is a non trivial problem to obtain the behavior of the Green 
function at s — 0. 

In the case of the Schwarzschild black hole, stability is demonstrated by 



Kay and Wald |117 who showed the boundedness of all solutions with data of 
compact support. 

The issue is more subtle for Kerr. There is a conserved energy, but because of 
the ergoregion its integrand is not positive definite, hence the conserved energy 
could be finite while the field still might grow exponentially in pa rts o f the 



spacetim e. P apers by Press and Teukolsky |165|, Hartle and Wilkins [109|, and 



Stewart |193| try to exclude the existence of an exponentially growing normal 



mode. Their work makes the stability very plausible but is not as conclusive as 
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the Wald-Kay result. However this is a delicate issue as we see if, for example, we 
multiply the Regge- Wheeler potential by a factor e: For any e > we obtain an 
infinite number of QNMs, for e = 0, however there is no QNM! Whiting [208| 
has proven that there are no exponentially growing modes, and in his proof 
he showed that the growth of t he modes is at most linear. Recent numerical 



evolution calculations [132, |l33t| for slowly and fast rotating Kerr black holes 
pick up all the expected features (QNM ringing, tails) and show no sign of 
exponential growth. It should be no ted t hat the massive scalar perturbations of 
Kerr are known to be unstable |TJ, |7^. These unstable modes are known 
to be very slowly growing (with growth times similar to the age of universe). 

Let us finally turn to the "completeness of QNMs" . A general mathematical 
theorem (spectral theorem) implies that for systems like strings or membranes 
the general solution can be expanded into a converging sum of normal modes. A 
similar result can not be expected for QNMs, the reasons are given in section ||. 
There is, however, the possibility that an infinite sum of the form (^5|) will be 
a representation of a solution for late times. This property has been shown 
by Beyer for the Poschl- Teller potential which has a similar form as the 
potential on Schwarzschild (23). The main difference is its exponential decay at 
both ends. In |158| NoUert and Price propose a definition of completeness and 
show its adcquateness for a particular model problem. There are also systematic 
studies about the relation between the structure of the QNM's of the Klein- 
Gordon equation and the form of the potential. In these studies there is a 
discussion on both the requirements for QNMs to form a complete set and the 
definition of completeness. 
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4 Quasi-Normal Modes of Relativistic Stars 



Pulsating stars are important sources of information for astrophysics. Nearly 
every star undergoes some kind of pulsation during its evolution from the early 
stages of formation until the very late stages, usually the catastrophic creation 
of a compact object (white dwarf, neutron star or black hole). Pulsations of 
supercompact objects are of great importance for relativistic astrophysics since 
these pulsations are accompanied by the emission of gravitational radiation. 
Neutron star oscillations were also proposed to explain the quasi-periodic vari- 



ability found in radio-pulsar and X-ray burster signals pOq , |146|] . In this chap- 
ter we shall discuss various features of neutron star non-radial pulsations i.e. 
the various modes of pulsation, mode excitation, detection probability and the 
possibility to extract information (to estimate, for example, the radius, mass 
and stellar equation of state) from the detection of the associated gravitational 
waves. It is not in our plans to discuss rotating relativistic stars; the interested 



reader should refer to another review in this journal [192|. Radial oscillations are 



also not discussed since they are not interesting for gravitational wave research. 

4.1 Stellar Pulsations: The Theoretical Minimum 

For the study of stellar oscillations we shall consider a spherically symmetric and 
static spacetime which can be described by the Schwarzschild solution outside 
the star, see equation (p7|). Inside the star, assuming that the stellar material 
is behaving like an ideal fluid, we define the energy momentum tensor 

T^u = {p + p)u^,Uv + PQtiu, (44) 

where p{r) is the pressure, p{r) is the total energy density. Then from the 
conservation of the energy-momentum and the condition for hydrostatic equi- 
librium we can derive the Tolman-Oppenheimer-Volkov (TOV) equations for 
the interior of a spherically symmetric star in equilibrium. Specifically, 

e-^ = l-^, (45) 
r 

and the "mass inside radius r" is represented by 

m(r) — An pr^dr. (46) 



Jo 

This means that the total mass of the star is M = ra{R), with R being the 
star's radius. To determine a stellar model we must solve 

rfp ^ P + P dv 

dr 2 dr' ^ ' 

where 

dv 2e^ {m + inpr^) 
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These equations should of course, be supplemented with an equation of state 
P = • • ■) as input. Usually is sufficient to use a one-parameter equation of 
state to model neutron stars, since the typical thermal energies are much smaller 
than the Fermi energy. The polytropic equation of state p = Kp^^'^/'^ where 
K is the polytropic constant and A'^ the polytropic exponent, is used in most of 
the studies. The existence of a unique global solution of the Einstein equations 
for a given equation of state and a given value of the central density has been 



proven by Rendall and Schmidt |]l74 | 



If we assume a small variation in the fluid or/and in the spacetime we must 
deal with the perturbed Einstein equations 

and the variation of the fluid equations of motion 

'^fe)=0, (50) 



while the perturbed metric will be given by equation (Jli 

Following the procedure of the previous section one can decompose the per- 
turbation equations into spherical harmonics. This decomposition leads to two 
classes of oscillations according to the parity of the harmonics (exactly as for 
the black hole case). The first ones called even (or spheroidal, or polar) produce 
spheroidal deformations on the fluid, while the second are the odd (or toroidal, 
or axial) which produce toroidal deformations. 

For the polar case one can use certain combinations of the metric pertur- 
bations as unknowns, and the linearized field equations inside the star will be 
equivalent to the following system of three wave equations for unknowns S, F, H: 



1 d'^S d'^S 



1 d'^F d^F 



1 d^H d^H 



+ Li{S,F,£)^0, (51) 
-L2{S,F,H,e) ^0, (52) 



(C3)2 dH d^r, 
and the constraint 



L3{H,H',S,S',F,F'J)^0, (53) 



d^F 



L4{F,F',S,S',H,e) ^0. (54) 



The linear functions Li, {i = 1, 2, 3, 4) depend on the background model and 
their explicit form can be found in | IIJ , |^ . The functions S and F correspond 



to the perturbations of the spacetime while the function H is proportional to 
the density perturbation and is only defined on the background star. With Cg 
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we define the speed of sound and witli a prime we denote differentiation witli 
respect to r*: 

^ = e(-^)/^^. (55) 
or* or 

Outside the star there are oniy perturbations of the spacetime. These are de- 
scribed by a single wave equation, the Z erilli equation mentioned in the previous 
section, see equations (|2l| ) and (^J). In |ll8[ it was shown that (for background 
stars whose boundary density is positive) the above system - together with the 
geometrical transition conditions at the boundary of the star and regularity 
conditions at the center - admits a well posed Cauchy problem. The constraint 
is preserved under the evolution. We see that two variables propagate along 
light characteristics and the density H propagates with the sound velocity of 
the background star. 



It is possible to eliminate the constraint - first done by Moncrief |152] - 
if one solves the constraint (|^ for H and puts the corresponding expression 
into ^2- (The characteristics for F change then to sound characteristics inside 
the star and light characteristics outside.) This way one has just to solve two 
coupled wave equations for S and F with unconstrained data, and to calculate 
H using the constraint from the solution of the two wave equations. Again the 
explicit form of the equation can be found in Q. 

Turning next to quasi-normal modes in the spirit of section ^ we can Laplace 
transform the two wave equations and obtain a system of ordinary differential 
equations which is of fourth order. The Green function can be constructed from 
solutions of the homogeneous equations (having the appropriate behavior at the 
center and infinity) and its analytic continuation may have poles defining the 
quasi-normal mode frequencies. 

From the form of the above equations one can easily see two limiting cases. 
Let us first assume that the gravitational field is very weak. The n eq uation ( ^T[ ) 
and (^2|) can be omitted (actually S* — + in the weak field limit ||200| , |j) and we 
find that one equation is enough to describe (with acceptable accuracy) the oscil- 
lations of the fluid. This approach is known as the Cowling approximation ||64| . 
Inversely, we can assume that the coupling between the two equations ( pl[ ) 
and ( |52|) describing the spacetime perturbations with the equation (^) is weak 
and consequently derive all the features of the spacetime perturbations from only 
the two of them. This is what is called the "inverse Cowling approximation" 
(ICA) 11. 

For the axial case the perturbations reduce to a single wave equation for the 
spacetime perturbations which describes toroidal deformations 

- + 5f + ^ W + 1)- + r\p-p)- 6M] = 0, (56) 

where X ~ hnj,. Outside the star, pressure and density are zero and this 
equation is reduced to the Regge- Wheeler equation, see equations ( pl| ) and (p^. 
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In Newtonian theory, if the star is non-rotating and the static model is a perfect 
fluid (i.e. shear stresses are absent), the axial oscillations are a trivial solution 
of zero frequency to the perturbation equations and the variations of pressure 
and density are zero. Nevertheless, the variation of the velocity field is not 
zero and produces non-oscillatory eddy motions. This means that there are no 



oscillatory velocity fields. In the relativistic case the picture is identical [ 202 | 
nevertheless; in this case there a re st ill QNMs, the ones that we will describe 
later as ^^spacetime or w-modes^^ [125|. 

When the star is set in slow rotation then the axial modes are no longer 
degenerate, but instead a new family of modes emerges, the so-called r-modes. 
An interesting property of these modes that has been pointed out by Anders- 
son Q is that these modes are generically unstable due the Chandrase khar- 



142| 



Friedman-Schutz instability |l53| , |95|| and furthermore it has been shown ]2 
that these modes can potentially restrict the rotation period of newly formed 
neutron stars and also that they can radiate away detectable amounts of grav- 
itational radiation [161|. The equations describing the perturbations of slowly 
rotating relativistic stars have been derived by Kojima |12C, 121 1, and Chan- 
drasekhar and Ferrari [ |6l| . 



4.2 Mode Analysis 

The study of stellar oscillations in a general relativistic context already has a 
history of 30 years. Nevertheless, recent results have shown remarkable features 
which had previously been overlooked. 

Until recently most studies treated the stellar oscillations in a nearly New- 
tonian manner, thus practically ignoring the dynamical properties of the space- 
time pol, |l9|, Im], lool |l^, Im^, ^, luel. The spacetime was used as 



the medium upon which the gravitational waves, produced by the oscillating 
star, propagate. In this way all the families of modes known from Newtonian 
theory were found for relativistic stars while in addition the damping times due 
to gravitational radiation were calculated. 

Inspired by a simple but instructive model |128], Kokkotas and Schutz 
showed the existence of a new family of modes: the w-modes [13C|. These are 
spacetime modes and their properties, although different, are closer to the black 
hole QNMs than to the standard fluid stellar modes. The main characteristics 
of the w-modes are high frequencies accompanied with very rapid damping. 
Furthermore, these modes hardly excite any fluid motion. The existence of these 



modes has been verified by su bsequent work 1 13S , 2J ; a part of the spectrum was 
found earlier by Kojima |ll£]and it has been shown that they exist also for odd 
parity (axial) oscillations [ |l25| . Moreover, sub- families of w-modes have been 
found for b oth t he polar and axial oscillations i.e. the interface modes found by 
Leins et al. [|l38[ (see also |p7|]), and t he trapped modes found by Chandrasekhar 



and Ferrari 1601] (see also |125, 122, |l^). Recently, it has been proven that 
one can reveal all the properties of the w-modes even if one "freezes" the fluid 
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mode 


frequency 


damping time 


/ 


2.87 kHz 


0.11 sec 


Pi 


6.57 kHz 


0.61 sec 


91 


19.85 Hz 


years 


Wl 


12.84 kHz 


0.024 ms 


WII 


8.79 kHz 


0.016 ms 



Table 2: Typical values of the frequencies and the damping times of various 
families of modes for a polytropic star (N — 1) with R — 8.86 km and M = 
1.27Mq are given, pi is the first p-mode, gi is the first g-mode ]^/ , wi stands 
for the first curvature mode and wn for the slowest damped interface mode. 
For this stellar model there are no trapped modes. 



oscillations (Inverse Cowling Approximation) ||2^. In the rest of this section 
we shall describe the features of both families of oscillation modes, fluid and 
spacetime, for the case £ = 2. 



4.2.1 Families of Fluid Modes 



For non-rotating stars the fluid modes exist only for polar oscillations. Here 
we will describe the properties of the most important modes for gravitational 
wave emission. These are the fundamental, the pressure and the gravity modes; 
this division has been done in a phenomenological way by Cowling ]6^ . For 
an extensive discussion of other families of fluid modes we refer the reader 
|99|| and [147, 146 1. Tables of frequencies and damping times of neutron 



to |g8|, 

star oscillations for twelve equations of state can be found in a recent work [[19| 
which verifies and extends earlier work |141|. In Table 2 we show characteristic 
frequencies and damping times of various QNM modes for a typical neutron 
star. 



• The /-mode (/mdamental) is a stable mode which exists only for non- 
radial oscillations. The frequency is proportional to the mean density of 
the star and it is nearly independent of the details of the stellar structure. 
An exact formula for the frequency can be derived for Newtonian uniform 
density stars 

This relation is approximately correct also for the relativistic case [0 
(see also the discussion in section 5.4). The /-mode eigenfunctions have no 
nodes inside the star, and they grow towards the surface. A typical neutron 
star has an /-mode with a frequency of 1.5 — 3 kHz and the damping time 
of this oscillation is less than a second (0.1 — 0.5 sec). Detailed data for 
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the frequencies and damping times (due to gravitational radiation) of the 



/-mode for various equations of state can be found in |141 , pi|. Estimates 
for the damping times due to viscosity can be found in ||69|, |7l[|. 

• The p-modes (pressure or acoustic) exist for both radial and non-radial 
oscillations. There are infinitely many of them. The pressure is the restor- 
ing force and it experiences substantial fluctuations when these modes are 
excited. Usually, the radial component of the fluid displacement vector is 
significantly larger than the tangential component. The oscillations are 
thus nearly radial. The frequencies depend on the travel time of an acous- 
tic wave across the star. For a neutron star the frequencies are typically 
higher than 4—7 kHz (pi-mode) and the damping times for the first few 
p-modes are of the order of a few seconds. Their frequencies and damping 
times increase with the order of the mode. Detailed data for the frequen- 
cies and damping times (due to gravitational radiation) of the pi-mode 
for various equations of state can be found in Q . 

• The 5-modes (gravity) arise because gravity tends to smooth out ma- 
terial inhomogeneities along equipotential level-surfaces and buoyancy is 
the restoring force. The changes in the pressure are very small along the 
star. Usually, the tangential components of the fluid displacement vec- 
tor are dominant in the fluid motion. The g-modes require a non-zero 
Schwarzschild discriminant in order to have non-zero frequency, and if 
they exist there are infinitely many of them. If the perturbation is stable 
to convection, the g- modes will be stable (w^ > 0); if unstable to con- 
vection the modes are unstable (w^ < 0); and if marginally stable to 
convection, the g-mode frequency vanishes. For typical neutron stars they 
have frequencies smaller than a hundred Hz (the frequency decreases with 
the order of the mode), and they usually damp out in time much longer 
than a few days or even years. For an extensive discussion about 5-modes 
in relativistic stars refer to e3 ' ^'^^ ^ study of the instability of the 



modes of rotating stars to gravitational radiation reaction refer to |134| 



The r- modes (rotational) in a non-rotating star are purely toroidal (axial) 
modes with vanishing frequency. In a rotating star, the displacement 
vector acquires spheroidal components and the frequency in the rotating 
frame, to first order in the rotational frequency of the star, becomes 

= WTTy ^''^ 

An inertial observer measures a frequency of 

LOi — ujj. — mfl. (59) 



From (58) and ( |59| ) it can be deduced that a counter-rotating (with re- 
spect to the star, as defined in the co-rotating frame) r-mode appears as co- 
rotating with the star to a distant inertial observer. Thus, all r-modes with 



26 



^ > 2 are generically unstable to the emission of gravitational radiation, 
due to the Chandrasekhar-Friedman-Schutz (CFS) mechanism [OSt - 
The instability is active as long as its growth-time is shorter than the 
damping-time due to the viscosity of neutron star matter. Its effect is to 
slow down, within a year, a rapidly rotating neutron star to slow rotation 
rates and this explains why only slowly rotating pulsars are associated with 
supernova remnants |23, 142| , 131 1. This suggests that the r-mode insta- 
bility might not allow millisecond pulsars to be formed after an accretion 
induced collapse of a white dwarf |2^. It seems that millisecond pulsars 
can only be formed by the accretion induced spin-up of old, cold neutron 
stars. It is also possible that the gravitational radiation emitted due to 
this instability by a newly formed neutron star could be detectable by the 
advanced versions of the gravitational wave detectors presently under con- 
struction [161|. Recently, Andersson, Kokkotas and Stergioulas [|^ have 
suggested that the r-instability might be responsible for stalling the neu- 
tron star spin- up in strongly accreting Low Mass X-ray Binaries (LMXBs). 
Additionally, they suggested that the gravitational waves from the neu- 
tron stars, in such LMXBs, rotating at the instability limit may well be 
detectable. This idea was also suggested by Bildsten and studied in 
detail by Levin |139[. 



4.2.2 Families of Spacetime or w-Modes 

The spectra of the three known families of w-modes are different but the spec- 
trum of each family is similar both for polar and axial stellar oscillations. As 
we have mentioned earlier they are clearly modes of the spacetime and from 
numerical calculations appear to be stable. 



The curvature modes are the standard w- modes |130|. They are the 
most important for astrophysical applications. They are clearly related 
to the spacetime curvature and exist for all relativistic stars. Their main 
characteristic is the rapid damping of the oscillations. The damping rate 
increases as the compactness of the star decreases: For nearly Newtonian 
stars (e.g. white dwarfs) these modes have not been calculated due to 
numerical instabilities in the various codes, but this case is of marginal 
importance due to the very fast damping that these modes will undergo. 
One of their main characteristics is the absence of significant fluid motion 
(this is a common feature for all families of w- modes). Numerical stud- 
ies have indicated the existence of an infinite number of modes; model 
problems suggest this too |12§| , |40| , [T^ . For a typical neutron star the 
frequency of the first w-mode is around 5—12 kHz and increases with the 
order of the mode. Meanwhile, the typical damping time is of the order 
of a few tenths of a millisecond and decreases slowly with the order of the 
mode. 



27 



4 

2 - 
10-' 
08 
06 - 
04 
02 - 



10- 
10-=T 
10-" 1 
10-S 
10- 



interface modes- 



0+0+0+' 



o+o+o 



curvature modes 



trapped modes 



O axial modes 
+ polar modes 
X black-hole 



50 75 1 00 

Re((BM) 



Figure 3: A graph which shows all the w-modes: curvature, trapped and 
interface both for axial and polar perturbations for a very compact uniform 
density star with MjR = 0.44. The black hole spectrum is also drawn for 
comparison. As the star becomes less compact the number o/ trapped modes 
decreases and for a typical neutron star (M/R = 0.2) they disappear. The 
Im{(jj) = 1/ damping of the curvature modes increases with decreasing com- 
pactness, and for a typical neutron star the first curvature mode nearly coincides 
with the fundamental black hole mode. The behavior of the interface modes 
changes slightly with the compactness. The similarity of the axial and polar 
spectra is apparent. 
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The trapped modes exist only for supercompact stars (i? < 3M) i.e. 
when the surface of the star is inside the peak of the gravitational field's 



potential barrier 125|. Practically, the first few curvature modes be- 



come trapped as the star becomes more and more compact, and even the 



/-mode shows similar behavior |122, [11. The trapped modes, as with 



all the spacetime modes, do not induce any significant fluid motions and 
there are only a finite number of them (usually less than seven or so). 
The number of trapped modes increases as the potential well becomes 
deeper, i.e. with increasing compactness of the star. Their damping is 
quite slow since the gravitational waves have to penetrate the potential 
barrier. Their frequencies can be of the order of a few hundred Hz to a 
few kHz, while their damping times can be of the order of a few tenths of 
a second. In general no realistic equations of state are known that would 
allow the formation of a sufficiently compact star for the trapped modes 
to be relevant. 



The interface modes |13S] are extremely rapidly damped modes. It 
seems that there is only a finite number of such modes (2 — 3 modes 
only) and they are in some ways similar to the modes for acoustic 
waves scattered off a hard sphere. They do not induce any significant fluid 
motion and their frequencies can be from 2 to 15 kHz for typical neutron 
stars while their damping times are of the order of less than a tenth of a 
millisecond. 



4.3 Stability 

The stability of radial oscillations for non-rotating stars in general relativity is 
well understood. Especially, the stability of static spherically symmetric stars 
can be determined by examining the mass-radius relation for a sequence of 



equilibrium stellar models, see for example Chapter 24 in |150|. The radial 
perturbations are described by a Sturm-Liouville second order equation with 
the frequency of the mode being the eigenvalue w^, then for real uj the modes 
will be stable while for imaginary w they will be unstable [M , see also Chapter 



17.2 in |18S 



The stability of the non-radially pulsating stars (Newtonian or relativistic) 
is determined by the Schwarzschild discriminant 

S[r) = ^l-^^J^, (60) 

where Fi is the star's adiabatic index. This can be understood if we deflne the 
local buoyancy force / per unit volume acting on a fluid element displaced a 
small radial distance 5r to be 

/ ^ -g{r)S{r)5r, (61) 
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where g is the local acceleration of gravity. When S is negative in some region 
the buoyancy force is positive and the star is unstable against convection, while 
when S is positive the buoyancy force is restoring and the star is stable against 
convection. Another way of discussing the stability is through the so-called 
Brunt- Vaisala frequency N"^ — gS{r) which is the characteristic frequency of 
the local fluid oscillations. Following earlier discussions when iV^ is positive, 
the fluid element undergoes oscillations, while when N'^ is negative the fluid is 
locally unstable. In other words, in Newtonian theory stability to non-radial 
oscillations can be guaranteed only if 5 > everywhere within the star . In 
general relativity Q, this is a sufficient condition, and so if 5 > the quasi- 
normal modes are stable. For an extensive discussion of stellar instabilities for 
both non-rotating and rotating stars (which are actually more interesting for 

the gravitational wave astronomy) refer to [177, 140, I92|. 

For completeness the same applies as outlined at the end of section 3.3. A 



model calculation of Price and Husain [168 , however indicated that the nearly 



Newtonian quasi-normal modes might be a basis for the fluid perturbations. 
Further mathematical investigation is needed to clarify this issue. 
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5 Excitation and Detection of QNMs 



A critical issue related to the discussions of the previous sections is the exci- 
tation of the QNMs. The truth is that, although the QNMs are predicted by 
our perturbation equations, it is not always clear which ones will be excited and 
under what initial conditions. As we have already mentioned in the introduc- 
tion there is an excellent agreement between results obtained from perturbation 
theory and full nonlinear evolutions of Einstein equations for head-on collisions 
of two black holes. Still there is a degree of arbitrariness in the definition of 
initial data for other types of stellar or black hole perturbations. This due to 
the arbitrariness in specifying the gravitational wave content in the initial data. 

The construction of acceptable initial data for the evolution of perturbation 
equations is not a trivial task. In order to specify astrophysically relevant ini- 
tial data one should first solve the fully nonlinear 3-dimensional initial value 
problem for (say) a newly formed neutron star that settles down after core col- 
lapse or two colliding black holes or neutron stars. Afterwards, starting from the 
Cauchy data on the initial hypersurface, one can evolve forward in time with the 
linear equations of perturbation theory instead of the full nonlinear equations. 
Then most of the long-time evolution problems of numerical relativity (throat 
stretching when black holes form, numerical instabilities or effects due to the 
approximate outer boundary conditions) are avoided. Additionally, the inter- 
pretations of the computed fields in terms of radiation is immediate This 
scheme has been used by Price and PuUin ]17C| | and Abrahams and Cook |Q 
with great success for head-on colliding black holes. The success was based on 
the fact that the bulk of the radiation is generated only in the very strong-field 
interactions around the time of horizon formation and the radiation generation 
in the early dynamics can be practically ignored (see also the discussion in |^ ) . 
The extension of this scheme to other cases, like neutron star collisions or su- 
pernovae collapse, is not trivial. But if one can define even numerical data on 
the initial hypersurface then the perturbation method will be probably enough 
or at least a very good test for the reliability of fully numerical evolutions. For a 
recent attempt towards applying the above techniques in colliding neutron stars 
see[||gl. 

Before going into details we would like to point out an important issue, 
namely the effect of the potential barrier on the QNMs of black holes. That is, 
for any set of initial data that one can impose, the QNMs will critically depend 
on the shape of the potential barrier, and this is the reason that the close limit 
approximation of the two black-hole collision used by Price and Pullin |17C| was 
so successful, because whatever initial data you provide inside the r < 3M region 
(the peak of the potential barrier is around 3A/) the barrier will "filter" them and 
an outside observer will observe only the QNM ringing (see for example recent 
studies by Allen, Camarda and Seidel 0). This point of view is complementary 
to the discussion earlier in this section, since roughly speaking even before the 
creation of the final black hole the common potential barrier has been created 
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and anything that was to be radiated had to be "filtered" by this common 
barrier. 

5.1 Studies of Black Hole QNM Excitation 

In the study of QNMs of black holes the attention, in most cases, was focused 
on estimating the spectrum and its properties. But there is limited work in 
the direction of understanding what details of the perturbation determine the 
strength of the QNM ringing. 

In 1977 Detweiler discussed the resonant oscillations of a rotating black 
hole, and after identifying the QNMs as "resonance peaks" in the emitted spec- 
trum he showed that the modes formally correspond to poles of a Green function 



to the inhomogeneous Teukolsky equation |197|. This idea has been extended 
in a more mathematically rigorous way by Leaver |136|. Leaver extracts the 
QNM contribution to the emitted radiation as a sum over residues. This sum 
arises when the inversion contour of the Laplace transform, which was used to 
separate the dependence on the spatial variables from the time dependence, is 
deformed analytically in the complex frequency plane. In t his way t he contri- 



bution from the QNM can be accounted for. Sun and Price [194, 195 1 discussed 



in detail the way that QNM are excited by given Cauchy data based to some 



extent on numerical results obtained by Leaver |136|. Lately, Andersson [y_2[ 
used the phase-integral method to determine some characteristics of the QNM 
excitation. 

5.2 Studies of Stellar QNM Excitation 

The discussion in the previous sections has shown that we have an acceptable 
knowledge of stellar pulsations and that we can extract information from the 
detection of such oscillations. But as was clear from the discussion earlier in 
this section, the energy released as gravitational radiation during the stellar 
collapse is basically unknown. Additionally, it is not known how the energy is 
distributed in the various fluid and spacetime modes. Both uncertainties depend 
strongly on the details of gravitational collapse, or in general the mechanism 
that excites the modes. Unless full 3D general relativistic codes are generated 
for the gravitational collapse or the final stages of binary coalescence we will 
never be able to give definite answers to the above questions. In the meantime, 
perturbation theory is a reliable way to get some first hints and indications. 
A survey of the literature reveals several indications that the pulsation modes 
are present in the gravitational wave signals from coalescing stars. Waveforms 



obtained by Nakamura and Oohara |155] show clear mode presence. Ruffert et 



al. |175| have also obtained gravitational wave signals from coalescing stars that 
show late-time oscillations. Their waveforms and spectra show oscillations at 
frequencies between 1.5 and 2 kHz. 
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Figure 4: The response of a Schwarzs child black hole as a Gaussian wave packet 
impinges upon it. The QNM signal dominates the signal after t ~ TOM while 
at later times (after t « 300MJ the signal is dominated by a power-law fall-off 
with time. 



The situation is similar for rotating core collapse. Most available studies use 
Newtonian hydrodynamics and account for gravitational wave emission through 
the quadrupole formula. The collapse of a non-rotating star is expected to 
bounce at nuclear densities, but if the star is rotating the collapse can also 
bounce at subnuclear densities because of the centrifugal force. In each case, 
the emerging gravitational waves are dominated by a burst associated with the 
bounce. But the waves that follow a centrifugal bounce can also show large 
amplitude oscillations that may be associated with pulsations in the collapsed 



core. Such results have been obtained by Monchmeyer et al. |15l| . Some of their 



models show the presence of modes with different angular dependence superim- 
posed. Typically, these oscillations have a period of a few ms and damp out in 
20 ms. The calculations also show that the energy in the higher multipoles is 
roughly three orders of magnitude smaller tha n that of the quadrupole. More 
recent simulations by Yamada and Sato |210| and Zwerger and Miiller [ 215 | 
also show post-bounce oscillations. The cited examples are encouraging, and it 
seems reasonable that besides the fluid modes the spacetime modes should also 
be excited in a generic case. To show that this is the case one must incorporate 
general relativity in the simulations of collapse and coalescence. As yet there 
have been few attempts to do this, but an interesting example is provided by the 



core collapse studies of Seidel et al. [187, 186, 182, 183 1. They considered axial 
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(odd parity) or polar (even parity) perturbations of a time-dependent back- 
ground (that evolves according to a specified collapse scenario). The extracted 
gravitational waves are dominated by a sharp burst associated with the bounce 
at nuclear densities. But there are also features that may be related to the fluid 
and the lu-modes. Especially for the axial case, since there are no axial fluid 
modes for a non-rotating star, it is plausible that this mode corresponds to one 
of the axial w-modes of the core. Furthermore, the power spectrum for one of 
the simulations discussed in |187| (cf. their Fig. 2) shows some enhancement 
around 7 kHz. 

Recently, Andersson and Kokkotas have studied the excitation of axial 
modes by sending gravitational wave pulses to hit the star (see figure ||). The 
results of this study are encouraging because they provide the first indication 
that the w-modes can be excited in a dynamical scenario. Similar results have 
recently been obtained by Borelli |Q for particles falling onto a neutron star. 
The excitation of the axial QNMs by test particles scattered by a neutron star 
have been recently studied in detail [203, 84, 27 and some interesting conclusions 
can be drawn. For example, that the degree of excitation depends on the details 
of the particle's orbit, or that in the cases that we have strong excitation of 
the quadrapole oscillations there is a comparable excitation of higher multipole 
modes. 

These results have prompted the study of an astrophysically more relevant 
problem, the excitation of even parity (or polar) stellar oscillations. In recent 
work Allen et al. |^ have studied the excitation of the polar modes using two 
sets of initial data. First, as previously for the axial case, they excited the 
modes by an incoming gravitational wave pulse. In the second case, the initial 
data described an initial deformation of both the fluid and the spacetime. In 
the flrst case the picture was similar to that of the axial modes i.e. the star was 
excited and emitted gravitational waves in both spacetime and the fluid modes. 
Nearly all of the energy was radiated away in the w-modes. 

This should be expected since the incoming pulse does not have enough time 
to couple to the fluid (which has a much lower speed of propagation of infor- 
mation). The infalling gravitational waves are simply affected by the spacetime 
curvature associated with the star, and the outgoing radiation contains an un- 
mistakable w-mode signature. In the second case one has considerable freedom 
in choosing the degree of initial excitation of the fluid and the spacetime. In the 
study the choice was some "plausible" initial data inspired by the treatment of 
the problem in the Newtonian limit. Then the energy emitted as gravitational 
waves was more evenly shared between the fluid and the spacetime modes, and 
one could see the /, p, and w modes in the signal. The characteristic signal was 
(as expected) a short burst (w-modes) followed by a slowly damped sinusoidal 
wave (fluid modes). 

The previous picture emerged as well in a recent study in the close limit of 
head-on collision of two neutron stars |^ . There the excitation of both fluid and 
spacetime modes is apparent but most of the energy is radiated via the fluid 
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Figure 5: Time evolution of axial perturbations of a neutron star, here only 
axial w-modes are excited. It is apparent in panel D that the late time behavior 
is dominated by a time tail. In the left panels (A and C) the star is ultra 
compact (M / R ~ Q.AA) and one can see not only the curvature modes hut also 
the trapped modes which damp out much slower. The stellar model for the 
panels (B and D) is a typical neutron star (M/R = 0.2) and we can see only 
the first curvature mode being excited. 
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modes. These new results show the importance of general relativity for the 
calculation of the energy emission in violent processes. Up to now the general 
way of calculating the energy emission has been the quadrupole formula. But 
this formula only accounts for the fluid deformations and motions, and as we 
have seen this accounts only for a part of the total energy. 




Figure 6: Excitation of polar modes due to an initial deformation of the star. 
The initial burst which is dominated mainly by the first w-mode is followed by a 
sinusoidal waveform dominated by the f and the first p-mode. In the upper panel 
the actual waveform is shotvn, while in the lower panel is its power spectrum. 
The wide peak in the power spectrum corresponds to the w-mode, while the sharp 
peaks correspond to the various fluid modes. 



5.3 Detection of the QNM Ringing 

It is well known in astrophysics that many stars will end their lives with a 
violent supernova explosion. This will leave behind a compact object which 
will oscillate violently in the first few seconds. Huge amounts of gravitational 
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radiation will be emitted and the initial oscillations will consequently damp out. 
The gravitational waves will carry away information about the compact object. 
If the supernova remnant is a black hole we will observe a short monochromatic 
burst lasting a few tenths of a millisecond (see figures [l]and|). If it is a neutron 
star we will observe a more complicated signal which will be the overlapping 
of several frequencies (see figure ||). The stellar signal will consist of a short 
burst (similar to that from a black hole) followed by a long lasting sinusoidal 
wave. Supernovae are q uite f requent, the expected event rate is 5.8 ± 2.4 events 
per century per galaxy [205|. The amount of energy emitted in such an event 
depends on the details of the collapse. A spherical collaps e wil l not produce any 
gravitational waves at all, while as much as lO^'^Af0C^ 11791 can be radiated 
away in a highly nonsphcrical one. The collapse releases an enormous amount of 
energy, at least equal to the binding energy of a neutron star, about O.ISMqc^. 
Most of this energy should be carried away by neutrinos, and this is supported 
by the neutrino observations at the time of supernova SN1987A. But even if 
only 1% of the energy released in neutrinos is radiated in gravitational waves 
then the above number makes sense. The present numerical codes used to 
simulate collapse predict that the energy emitted as gravitational waves will 
be of the order of 10"'* — 1O~^M0C^ However, most of these codes are 

based on Newtonian dynamics and the few fully general relativistic ones are 
not 3-dimensional. Modern computers are still not able to perform realistic 
simulations of gravitational collapse in 3D, including all the important nuclear 
reactions and neutrino and photon transport. For example, most of the codes 
fail to explain the high average pulsar velocity which is believed to be a result 
of a boost that the neutron star gets during the collapse due to anisotropy in 
the neutrino distribution | pT| . 

Although collapse may be the most frequent source for excitation of black 
hole and stellar oscillations there are other situations in which significant pulsa- 
tions take place. For example, after the merger of two coalescing black holes or 
neutron stars it is natural to expect that the final object will oscillate. Thus the 
well known waveform for inspiralling binaries |p5| will be followed by a short, 
but not yet properly known, period (the merger phase) and will end with the 
characteristic quasi-normal signal (ringing) of the newly created neutron star or 
black hole. During the inspiralling pha se the stellar oscillations can be excited 
by the tidal fields of the two stars |127| . A detailed description of the gravita- 
tional wave emission and detection from binary black hole coalescences can be 
found in two recent articles by Flanagan and Hughes ||9^, |^ . In the same way 
smaller bodies falling on a neutron star or black hole will excite oscillations. 
Stellar or black hole osci llati ons can also be excited by a close encounter with 
another compact object |203, |2^ . 

Another potential excitation mechanism for stellar pulsation is a starquake, 
e.g., associated with a pulsar glitch. The typical energy released in this process 
may be of the order of 10^^'^ Mqc^ . This is an interesting possibility considering 
the recent discovery of so-called magnetars: Neutron stars with extreme mag- 
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netic fields These objects are sometimes seen as soft gamma-ray repeaters, 
and it lias been suggested that the observed gamma rays are associated with 
starquakcs. If this is the fraction of the total energy could be released 

through nonradial oscillations in the star. As a consequence, a burst from a soft 
gamma-ray repeater may be associated with a gravitational wave signal. 

Finally, a phase-transition could lead to a sudden contraction during which 
a considerable part of the stars gravitational binding energy would be released, 
and it seems inevitable that part of this energy would be channeled into pulsa- 
tions of the remnant. Transformation of a neutron star into a strange star is 
likely to induce pulsations in a similar fashion. 

One way of calibrating the sensitivity of detectors is to calculate the ampli- 
tude of the gravitational wave that would be produced if a certain fraction of 
the released energy were converted into gravitational waves. To obtain rough 
estimates for the typical gravitational wave amplitudes from a pulsating star we 
use the standard relation for the gravitational wave flux which is valid far away 
from the star ]17^ 

where h is the gravitational wave amplitude and r the distance of the detector 
from the source. Combining this with i) dE/dt = E/2t where r is the damping 
time of the pulsation and E is the available energy, ii) the assumption that 
the signal is monochromatic (with frequency /) , and iii) the knowledge that the 
effective amplitude achievable after matched filtering scales as the square root of 
the number of observed cycles, icff = h^/n = h^/Jr, we get the estimates |18, 19| 

for the /-mode, and 



eff ^ 9.7 X 10-22 (^^_^) irr^l (64) 



E y^VlO kHzy^VsO kpc 

for the fundamental w-mode. Here we have used typical parameters for the 
pulsation modes, and the distance scale used is that to SN1987A. In this volume 
of space one would not expect to see more than one event per ten years or so. 
However, the assumption that the energy release in gravitational waves in a 
supernova is of the order of 10~^Mqc^ is very co nserv ative [179|. 
Similar relations can be found for black holes \l7i]: 

for stellar black holes, and 
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for galactic black holes. 

An important factor for the detection of gravitational waves are the pulsation 
mode frequencies. Existing resonant gravitational wave detectors, as well as 
laser interferometric ones which are under construction, are only sensitive in a 
certain bandwidth. The spherical and bar detectors are typically tuned to 0.6 — 
3 kHz, while the interferometers are sensitive within 10 — 2000 Hz. The initial 
part of the QNM waveform, which carries away whatever deformation a collapse 
left in the spacetime, is expected to be for a neutron star in the frequency range 
of 5 — 12 kHz (w-mode). The subsequent part of the waveform is constructed 
from combination of the /- and p-modes. Still the present gravitational wave 
detectors are sensitive only in the frequencies of the /-mode. For a black hole 
the frequency will depend on the mass and rotation rate^, thus for a 10 solar 
mass black hole the frequency of the signal will be around 1 kHz, around 100 Hz 
for a 100 Mq black hole and around 1 mHz for galactic black holes. 



5.4 Parameter Estimation 



For astronomy it is important not only to observe various astronomical phe- 
nomena but also to try to mine information from these observations. From 
the observations of solar and stellar oscillations (of normal stars) astronomers 
have managed to get details of the internal structure of stars. In our days the 
GONG program [110| for detailed observation of the solar seismology is well 
underway. This has suggested that, in a similar way, information about neu- 
tron star parameters (mass, radius), and internal structure or the mass and the 
rotation rate of black holes can be found, using the theory of QNMs. It will 
be instructive to briefly examine the case of oscillating black holes since they 
are much "cleaner" objects than stars. From the normal mode analysis of black 
hole oscillations we can get a spectrum which is related to the parameters of the 
black hole (mass M and angular momentum a). In particular, for the frequency 
of the first quasi-normal mode (which as we have stated previously is the most 
important one for the gravitational wave detection) the following approximate 
relations have been suggested |83, p9| : 
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These two relations can be inverted and thus from the "observed" frequency 
and the damping time we can derive the parameters of the oscillating black 
hole. In practice, the noise of the detector will contaminate the signal but still 

^For neutron stars the frequencies depend not only on the mass and rotation rate, but also 
on the radius and the equation of state. 
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(depending on the signal to noise ratio) we will get a very accurate estimate of 
the black hole parameters. A similar set of empirical relations cannot be derived 
in the case of neutron star oscillations since the stars are not as "clean" as black 
holes, since more than one frequency contributes. Although we expect that most 
of the dynamical energy stored in the fluid oscillations will be radiated away 
in the /-mode, some of the p-modes may be excited as well and a significant 
amount of energy could be radiated away through these modes (|] . As far as the 
spacetime modes are concerned we expect that only the curvature modes (the 
standard w-modes) will be excited, but it is possible that the radiated energy 
can be shared between the first two w-modes Nevertheless, Andersson and 
Kokkotas using the properties of the various families of modes (/, p, and 
w), managed to create a series of empirical relations which can provide quite 
accurate estimates of the mass, radius and equation of state of the oscillating 
star, if the / and the first w-mode can be observed. In figure Q one can see 
an example of the relation between the stellar parameters and the frequencies 
of the / and the first w-mode for various equations of state and various stellar 
models. There it is apparent that the relation between the /-mode frequencies 
and the mean density is almost linear, and a linear fitting leads to the following 
simple relation: 



w/(kHz) « 0.78 1.635 
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(69) 



We can also find the following relation for the frequency of the first w-mode: 
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(70) 



From tests performed using polytropic stars to provide data for the above 
relations it was seen that these equations predict the masses and the radii of 
the polytropes usually with an error less than 10%. There is work underway 
towards extracting the parameters of the star from a noisy signal |12(:]. 
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Figure 7: The upper graph shows the numerically obtained /-mode frequencies 
plotted as functions of the mean stellar density. In the second graph the func- 
tional RlUw is plotted as a function of the compactness of the star {M and R 
are in km, w/.modc and uJw-modc in kHz). The letters A, B, C, . . . correspond to 
different equations of state for which one can refer to [M . 
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6 Numerical Techniques 



For each candidate source of gravitational waves, gravitational wave astronomy 
needs answers to two questions. Firstly, how much energy will be carried away by 
the emitted gravitational waves? Secondly, which are the "preferred" frequencies 
at which a 10 Mq black hole or a neutron star will oscillate? The answer to 
the first question is that the energy will depend on the degree of asymmetry 
that the process generates, and it will depend critically on the initial data. 
In the previous section we have tried to provide some guesstimates for the 
energy emitted during the oscillation phase of black holes and neutron stars. 
The answer to the second question is related to the numerical solution of the 
perturbation equations. The numerical schemes developed for this purpose will 
be described in this section. 

Let us describe why the numerical calculation of quasi-normal mode fre- 
quencies is delicate. Consider again the case treated in section || of the wave 
equation with a potential with compact support. We try to find a complex 
number s with negative real part such that the solution which is e^'*^ for large 
positive a;, is e*^ for large negative x. Note that these solutions grow expo- 
nentially with \x\ and therefore one has to be very careful to make sure that 
there is no exponentially decaying part in the solution. The situation becomes 
even more complicated if we do not know explicitly because one can not 
characterize the correct solution by some growth property. This is for example 
the case for the Schwarzschild solution. 

6.1 Black Holes 

We would like to point out here that the attempts to calculate the QNM fre- 
quencies date back to the beginning of 70s. More specifically in their study of 
the black hole oscillations excited by an infalling particle Davis et al. [|7^ found 
that the peak of the spectrum is (for a Schwarzschild black hole) at around 
Mu) = 0.32 (geometrical units). This number is very close to the one calculated 
by much more accurate methods later. In what will follow we will describe the 
various methods used to calculate the QNM frequencies. 

6.1.1 Evolving the Time Dependent Wave Equation 

This approach was actually the first one used to study the QNM excitation by 
Vishveshwara [ p07[ but it has only been recently revived thanks to increased 
power of computers. In general, one does not need to Fourier decompose the 
perturbed Einstein equations but instead evolve them for given sets of initial 
data and at the end to Fourier analyze the resulting waveform. This procedure 
has certain advantages since one does not need to be so careful in considering 
the appropriate boundary conditions on the horizon, at infinity, on the surface 
or at the center of the star. This does not mean that these boundaries are 
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not important for the evolution schemes. The difference is that in the time 
independent case one formulates a boundary value problem, and the eigenvalues 
(quasi-normal frequencies) depend critically on the correct conditions on the 
various boundaries. A major disadvantage of the evolution schemes is that one 
cannot get the complete spectrum of the QNMs neither for a star nor for a black 
hole. The reason is that although any perturbation is the sum of the harmonics 
involved, in practice only a few of them will be observed and in the best case one 
can succeed in getting a few extra modes by "playing" around with the initial 
data. 

To be more specific, by evolving a perturbation on a black hole background 
one can get a QNM signal as the one shown in figure ^, but in this signal the 
Fourier transform will show that there are present at most two frequencies (the 
slowest damped ones); then by doing "matched filtering" of this signal we will 
get the right frequencies and damping times, but then all the extra modes of 
the spectrum are missing! See for example the work by Bachelot and Motet- 



Bachelot and Krivan et al |132, 132]. This is true also for stars; in 



recent evolutions of axial perturbations of stellar backgrounds Andersson and 
Kokkotas ||l8| saw only a few of the w-modes (the ones that damped slowest), 
while in similar calculations for even parity stellar perturbations Allen et al. |^ 
have seen only the /-mode, 2 — 3 p-modes and two of the w-modes. 

Of course with more detailed studies for various sets of initial data one might 
be successful to get a few more modes, but more important than deriving extra 
modes is understanding the physical situation which generates the appropriate 
initial data. 

Finally, the evolutions of the time dependent perturbation equations can 
be extremely useful (and probably will be the only way) for the calculation of 
the QNM frequencies and waveforms for the perturbations of the Kerr-Newman 
black hole and for slowly and fast rotating relativistic stars. 

6.1.2 Integration of the Time Independent Wave Equation 

This technique was used by Chandrasekhar and Detweiler |5|] and is based on 
the definition of QNMs given in section ^. They assumed that a QNM is a so- 
lution corresponding to incoming waves on the horizon and outgoing at infinity. 
Then by taking a series expansion of the Zerilli wave equation (^ |2^, ^) at 
both limits (horizon and infinity) of the form given by ( p7| ) they found initial 
values for the numerical integration of the equation. Their integration goes 
from both limits towards a common point which was set close to the peak of 
the potential i.e. around r = 3M. The values of cu for which the Wronskian of 
the two numerically taken solutions vanishes are the quasieigensolutions of the 
problem. In this way they managed to calculate the first 2 — 3 QNM frequencies 
of the Schwarzschild black hole for various harmonic indices. The accuracy of 



the method improves for increasing £. Later, Gunter 108 and Kokkotas with 



Schutz 129 used the same approach to calculate the QNMs of the Reissner- 
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Nordstrom black hole. 

The approach used by Nollert and Schmidt |156, 159| is more elaborate and 
based on a better estimate of the values of the quasi-eigenfunctions on both 
boundaries (±00); this leads to a more accurate estimate of frequencies and one 
also finds frequencies which damp extremely fast. Andersson |^ suggested an 
alternative integration scheme. The key idea is to separate ingoing and outgoing 
wave solutions by numerically calculating their analytic continuations to a place 
in the complex r-coordinate plane where they have comparable amplitudes. This 
method is extremely accurate. 



6.1.3 WKB Methods 



This technique, originally in the form suggested by Schutz and Will [180|, based 
on elementary quantum mechanical arguments, was later developed into a pow- 
erful technique with which accurate results have been derived. The idea is that 
one can reduce the QNM problem into the standard WKB treatment of scatter- 
ing of waves on the peak of the potential barrier. The simplest way to find the 
QNM frequencies is to use the well known Bohr-Sommerfeld (BS) rule. Using 
this rule it is possible to reproduce not only the Schutz- Will formula ( |30| ) but 
also to give a wa y to extend the accuracy of that formula by taking higher order 



terms 



114 , |105|| . The classical form of the BS rule for equations like ( pT| ) is 



- V{r)\ 



1/2 



dr = (n + 2^''^' 



(71) 



where rA,rB are the two roots (turning points) of oj^ — V{r) = 0. A more 
general form can be found ^% |4l} | which is valid for complex potentials. This 
form can be extended to the complex r plane where the contour encircles the 
two turning points which are connected by a branch cut. In this way one can 
calculate the eigenfrequencies even in the case of complex potentials, as it is the 
case for Kerr black holes. 

This method has been used for the calculation of the eigenfrequencies of 



the Schwarzschild [113|, Reissner-Nordstrom [129|, Kerr [185, 123| and Kerr- 
Newman |124| black holes (restricted case). In general with this approach one 
can calculate quite accurately the low-lying (relatively small imaginary part) 
QNM modes, but it fails to give accurate results for higher-order modes. 

This WKB approach was improved considerably when the phase integral 
formalism of Froman and Froman |]97|] was introduced. In a series of papers [p6| , 
the method was developed and a great number of even extremely 
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fast damped QNMs of the Schwarzschild black hole have been calculated with 
a remarkable accuracy. The application of the method for the calculation of 
QNM frequencies of the Reis sner-Nordstrom black hole [|6) has considerably 
improved earlier results |l29|] for the QNMs which damp very fast. Close to 
the logic of this WKB approach were the attempts of Blomc, Mashhoon and 
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Ferrari |4^, Q to calculate the QNM modes using an inversion of the black 
hole potential. Their method was not very accurate but stimulated future work 
using semianalytic methods for estimating QNMs. 

6.1.4 The Method of Continued Fractions 



In 1985, Leaver |135| presented a very accurate method for calculating the QNM 
frequencies. His method can be applied to the calculation of the QNMs of 
Schwarzschild, Kerr and, with some modifications, Reissner-Nordstrom black 



holes 1 137 . This method is very accurate also for the high-order modes. 

His approach was analogous to the determination of the eigenvalues of the H+ 
ion developed in [Q. A series representation of the solution /~ is assumed to 
represent also for the value of the quasi-normal mode frequency. For normal 
modes the method may work because is certainly bounded at infinity. In the 
case of quasi-normal modes this is not so clear because /"*" grows exponentially. 
Nevertheless, the method works very well numerically and was improved by 



NoUert [157| such that he was able to calculate very high mode numbers (up 
to 100,000!). In this way he obtained the asymptotic distributions of modes 
described in (|3l|). An alternative way of using the recurrence relations was 



suggested i n |145 | 



Nollert [|156[ explains in his PhD thesis why the method works. As initiated 



by Heisenberg et al. |111| he considers potentials depending analytically on a 
parameter A such that for A = 1 the potential has bound states - normal modes 
- and for for A = — 1 just quasi-normal modes. This is, for example, the case if 
we multiply the Regge- Wheeler potential ( p3| ) by — A. Assuming that the modes 
depend continuously on A, one can try to relate normal modes to quasi-normal 
modes and their methods of calculation. 

In the case of QNMs of the Kerr black hole, one has to deal in practice with 
two coupled equations, one which governs the radial part ( ^0| ) and another which 
governs the angular dependence of the perturbation (^9|). For both of them one 
can construct recurrence relations for the coefficients of the series expansion of 
their solutions, and through them calculate the QNM frequencies. 

For the case of the QNMs of the Reissner-Nordstrom black hole, the asymp- 
totic form of the solutions is similar to that shown in equation ( |28| ) but the 
coefficients a„ are determined via a four term recurrence relation. This means 
that the nice properties of convergence of the three term recurrence relations 
have been lost and one should treat the problem with great caution. Neverthe- 



less, Leaver [137 has overcome this problem and showed how to calculate the 
QNMs for this case. 

As a final comment on this excellent method we should point out that it 
has a disadvantage compared to the WKB based methods in that it is a purely 
numerical method and it cannot provide much intuition about the properties of 
the QNM spectrum. 
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6.2 Relativistic Stars 



Although the perturbation equations in the exterior of a star are similar to 
those of the black-hole and the techniques described earlier can be applied here 
as well, special attention must be given to the interior of the star where the 
perturbation equations are more complicated. 

For the time independent case the system of equations p^ , |3|) inside 
the star reduces to a 4th order system of ODEs |Q . One can then even treat 
it as two coupled time independent wave equations^. The first equation will 
correspond to the fluid and the second equation will correspond to spacetime 
perturbations. In this way one can easily work in the Cowling approximation 
(ignore the spacetime perturbations) if the aim is the calculation of the QNM 
frequencies of the fluid modes (/, p, g, . . .) or the Inverse Cowling Approxi- 
mation (ignore the fluid perturbations) if the interest is in w-modes. The 
integration procedure inside the star is similar to those used for Newtonian stars 
and involves numerical integration of the equations from the center towards the 
surface in such a way that the perturbation functions are regular at the center 
of the star and the Lagrangian variation of the pressure is zero on the surface 



(for more details refer to 141 , 13C]). The integrations inside the star should 
provide the values of the perturbation functions on the surface of the star where 
one has to match them with the perturbations of the spacetime described by 
Zerilh's equation (|l|, ||, ||). 

In principle the integrations of the wave equation outside the star can be 
treated as in the case of the black holes. Leaver's method of continued fraction 



has been used in 1 119113? | , Andersson's technique of integration on the complex 



r plane was used in l2l[] while a simple but effective WKB approach was used 



by Kokkotas and Schutz ||l30|, |211 | 



Finally, there are a number of additional approaches used in the past which 
improved our understanding of stellar oscillations in GR. In the following para- 
graphs they will be discussed briefly. 



Resonance Approach. This method was developed by Thorne |199|, 
the basic assumption being that there are no incoming or outgoing waves 
at infinity, but instead standing waves. Then by searching for resonances 
one can identify the QNM frequencies. The damping times can be esti- 
mated from the half-width of each resonance. This is a simple method 
and can be used for the calculations of the fluid QNMs. In a similar fash- 
ion Chandrasekhar and Ferrari ||5^ have calculated the QNM frequencies 
from the poles of the ratio of the amplitudes of the ingoing and outgoing 



* Chandrasekhar and Ferrari ]59| have also reduced the time independent perturbation 
equations (using a different gaugeyinto a 5th order system which involves only the spacetime 
perturbations with the fluid perturbations being calculated via a lgebraic r elations from the 



spacetime perturbations. It was later proven by Ipser and Price |112, 169| that this system 
of ODEs can be reduced to the standard eq uations in the Regge- W heeler gauge. That this is 
possible is apparent from equations |52| , p3| and it is discussed in |Q 
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waves. 



Direct Numerical Integration. This method was used by Lindblom 



and Detweiler |141] for the calculation of the frequencies and damping 
times of the /-modes for various stellar models for thirteen different equa- 
tions of state. In this case, after integration of the perturbation equations 
inside the star, one gets initial data for the integration of the Zerilli equa- 
tion outside. The numerical integration is extended up to "infinity" (i.e. 
at a distance where the solutions of Zerilli equations become approxi- 
mately simple sinusoidal ingoing and outgoing waves), where this solution 
is matched with the asymptotic solutions of the Zerilli equations which 
describe ingoing and outgoing waves. The QNM frequencies are the ones 
for which the amplitude of the incoming waves is zero. This method is 
more accurate than the previous one at least in the calculation of the 
damping times as has been verified in , but still is not appropriate for 
the calculation of the w-modes. 

Variational Principle Approach. Detweiler and Ipser derived a 
variational principle for non-radial pulsational modes. Associated with 
that variational principle is a conservation law for the pulsational energy 
in the star. The time rate of change of that pulsational energy, as given 
by the variational principle, is equal to minus the power carried off by 
gravitational waves. This method was used widely for calculating the 
/ Q and .g-modes ||^ and in studies of stability ||^, |7|] . 

WKB. This is a very simple method but quite accurate, and contrary 
to the previous three methods it can be used for the calculation of the 
QNM frequencies of the w-modes (this was the first method used for the 
derivation of these modes). In practice one substitutes the numerical 
solutions of the Zerilli equation with their approximate WKB solutions 
and identifies the QNM frequencies as the values of the frequency for 



which the amplitude of the incoming waves is zero [130, 211 



WKB-Numerical. This is a combination of direct integration of the 
Zerilli equation and the WKB method. The trick is that instead of in- 
tegrating outwards, one integrates inwards (using initial data at infinity 
for the incoming wave solution); this procedure is numerically more stable 
than the outward integration. On the stellar surface one needs of course 
not the solution for incoming waves but the one for outgoing ones. But 
through WKB one can derive approximately the value of the outgoing 
wave solution from the value of the incoming wave solution. In this way 



errors are introduced, nevertheless the results are quite accurate [130 
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7 Where Are We Going? 



From the previous discussions we can draw some general conclusions and sug- 
gestions for future work. We believe that independent of the advancement of 
numerical relativity and computer power there remains much future work in 
perturbation theory, in particular parallel to the numerical relativity efforts. 



7.1 Synergism Between Perturbation Theory and Numer- 
ical Relativity 

Fifteen years ago with the advancements in computer power, we could not think 
that perturbation theory would continue to play such an important role in many 
problems. Today it is widely accepted that there is a need for the development 
of new approximation schemes to accompany large scale simulations. Fully 
relativistic computer simulations give only numerical answers to problems; of- 
ten these answers do not provide physical understanding of which principles 
are important, or even what principles govern a given process. Even more, in 
some cases, simulation results can be simply incorrect or misleading. By closely 
coupling various perturbation schemes it is possible to interpret and confirm 
simulation results. 

For example, during the late stages of black hole or stellar coalescence or 
supernovae collapse, the system settles down to a slightly perturbed black hole 
or neutron star. Numerical codes, evolving the full nonlinear Einstein equations, 
should be able to accurately compute the waveforms required for gravitational 
wave detection. Parallel to this, it should be possible to evolve the perturbations 
of both black holes and stars (governed by their own linear evolution equations) 
for the same set of initial data. This is an important check of the fully nonlinear 
codes. An excellent example is the head-on collision of two black holes. This 
sounds like an impossible task for perturbation theory but it can be achieved 
if the two black holes are close together and they can be considered as having 
already merged into a single perturbed black hole (close limit); see more details 



in a recent review by PuUin |172]. Much work has already been performed for 
head-on collisions of two non-rotating black holes but the more realistic problem 
for the inspiral collision of rotating black holes is still open. 

Following in the same spirit are more recent calculations |^ of the close limit 
for two identical neutron stars. The results show the excitation of stellar QNMs 
and it remains for numerical relativity to verify the results. 

Finally, perturbation theory can be also used as a tool to construct a gauge 
invariant measure of the gravitational radiation in a numerically generated per- 
turbed black hole spacetime |@, 184 1. 



48 



7.2 Second Order Perturbations 



A basic feature of linearized perturbation theory is that there is no "built-in" 
indication of how good the approximation is. But if one has results for a phys- 
ical quantity to second order in a perturbation parameter, then the difference 
between the results of the first and second order theory is a quantitative in- 
dication of the error in the perturbation calculation. From this point of view 
second order perturbation theory is a practical tool for calculations. Neverthe- 
less, there remain many technical problems in this approach, for example the 

50[ | and the amount of calculations. Thus second 
order perturbation theory has turned out to be much more difficult than lin- 
earized theory, but if one overcomes these difhculties second order calculations 
will be a great deal easier than numerical relativity. From this point of view 
we encourage work in this direction. In particular, the perturbations of stellar 
oscillations should be extended to second order, and the study of the second 
order perturbations of black holes extended to the Kerr case. 



gauge that should be used |10C, 



7.3 Mode Calculations 

Although the modes have been well studied there remain a few open issues. In 
the black hole case, the Kerr-Newman spectrum is known only in a restricted 



case [124], so the general case should be studied, probably by evolving the 
perturbation equations. For non-rotating stars, there remain certain technical 
questions, i.e. to find if there exist a class of w-modes with even larger imaginary 
part or the existence of extra interface modes. But for rotating stars, even slowly 
rotating ones, there remains much work |192( |. There should be estimates of the 
fluid and spacetime modes for slowly rotating stars for various realistic equations 
of state, and the results should be incorporated into the method suggested in 
for the estimates of the stellar parameters. There should be work towards 



understanding the possible interaction of r- modes with 5- modes | 176 , |93| . And 
finally the time dependent perturbation equations for rotating stars should be 
evolved, because in this way we expect to see most of the new features that 
rotation induces in the spectra (splitting, instabilities etc). 



7.4 The Detectors 

As we have mentioned earlier in section ||, there are techniques available for 
the extraction of the QNM signal from the noise of the detectors |9^, ^ . 
But there are still issues related to the sensitivity of the planned detectors in 
parts of the spectrum. For example, the QNM frequencies of stellar black holes 
(10— IOOMq) will be of around 100— 1000 Hz, i.e. in the frequencies where the 
laser interferometers are sensitive. The QNM frequencies of galactic size black 
holes will be detectable only from space (LISA) since their frequencies will be in 
the mHz regime. For the QNM frequencies of stars, there is a lack of appropriate 
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detectors. Laser interferometers will be sensitive enough in the frequency regime 
of the /-mode but it will be very hard to detect signals in the frequencies of the 
p- and ly-modes. Nevertheless, from the discussion in section ^ it is apparent 
that there is a wealth of information in the signal of oscillating neutron stars, 
and in order to extract this information we need the p- or/and w-modes. This 
suggests that the ideas considering detectors, or arrays of detectors, in this high 
frequency regime should be considered more seriously. 
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9 Appendix: Schrodinger Equation Versus Wave 
Equation 

The purpose of this appendix is a brief comparison of properties of the Schro- 
dinger equation 

+ (72) 



and the wave equation 



2 



$=(-^+y(x))$ (73) 



for the same potential V{x). 

The ansatz of "stationary states" 

^^J = e-'^'i^{x) (74) 

leads to the time independent Schrodinger equation 

- V'" + V" = E^j. (75) 

The ansatz (corresponding to Laplace transformation) 

$ = e''(t){x) (76) 

gives for the wave equation 

-(j)" + V (t) = -s^(t). (77) 



The equations (|7^) and (77) are the same if we set E — ~s . Hence time 
independent scattering theory - the theory of the operator ~cP /dx"^ + V{x) 
— does not know whether we deal with the wave equation or the Schrodinger 
equation. 

Consider first a positive potential of compact support. From quantum me- 
chanics we know that there are no bound states. This is easy to understand: 
outside the support of the potential the solutions are exp{zLx^/ —E), hence only 



negative values of E are possible eigenvalues. However, if we multiply (75) hy tJj 
and integrate over all x we obtain a contradiction because of the positivity of V . 
Hence the operator —(P/dx"^ -\- V{x) has no eigenfunctions and there are only 
scattering states, the continuous spectrum. The operator is selfadjoint on the 
Hilbert space of square integrable functions on the real line. Its resolvent, R^; 
is defined for all complex E outside the continuous spectrum, which consists of 
the non negative real numbers. R^; is an integral operator whose kernel is the 
Green function constructed in section^. There G(s,x,x') was considered as a 
function of s. As a function of -E, G is analytic on the whole complex plane with 
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the exception of real E < 0, the continuous spectrum. In terms of s, the Green 
function - and the resolvent - is defined on the "physical half space" Re{s) > 0. 

Resonances of Schrodinger operators can be defined as the poles of the ana- 
lytic extension of the Green function or the resolvent. There is a huge amount 
of literature on the subject, in particular mathematical papers. A convenient 
starting point may be the proceedings of a conference on resonances in 1984 Q. 
Existence of resonances, asymptotic distribution and also their interpretation is 
treated. There is in particular the recent development of "Geometric scattering 
Theory" , where the use of "pseudo differential operators" gives strong and inter- 



esting results [ |l49[ . We describe one such result on the asymptotic distribution 
of quasi-normal mode frequencies of the Schwarzschild spacetime in section ||. 
It is amusing to note that in the field of quantum mechanics the same difficulties 
in defining the notion of a resonance occurred as in relativity in the context of 



"normal modes of black holes"! |189| is a good reference explaining this point. 

Historically, quasi-normal modes appeared the first time in Gamow's treat- 
ment of the Qf-decay. The model he studies is a potential with two positive 
square potentials. Radioactive decay is exponential in experiments. However, 
even the decay in time of solutions of the free Schrodinger equation is a power 
law decay (l/t for 1-dimensional systems), similarly for potentials with com- 
pact support. So we face the difficulty to characterize that part of the time 
evolution, in which exponential decay is a good approximation, knowing that 
the final decay is polynomial! In ]19C| | such estimates are derived. 

Let us finally consider general potentials of compact support. If the poten- 
tial well is deep enough a finite number of negative eigenvalues £"„ may exist 
describing bound states of the quantum system. What are the properties of the 
corresponding solutions of the wave equation? Let ipni^) be an eigenfunction 
of -(f /dx^ + y{x) with eigenvalue i?„ < 0. Then (E = -s^) 

$ = e^^Vn(:c) (78) 

is a solution of the wave equation. For large positive x we have 

= e-^^ ^ (79) 

as the only square integrable solution of the Schrodinger equation with vanishing 
potential. Thus the solution of the wave equation for large x is 



-E„ (t-x) 



(80) 



This solution grows exponentially in time and falls off exponentially in space 
for X — > oo. This apparently strange behaviour is possible because the energy 
density of the conserved energy of the wave equation 

e= ^ + + ^*') (81) 
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is not positive definite, if the potential is somewhere negative. In this situation 
the solution can grow in time and nevertheless have conserved finite energy. 

In the complex s plane the eigenvalues appear as poles of the Green func- 
tion with values s„ = ^—En- To the right of the largest eigenvalue we have 
analyticity of the Green function. 

Let us close this section by remarking that functional analysis techniques 
can also be used to develop existence theory: In the case discussed above, 
—(P/dx^ + V{x) is selfadjoint on the space of square-integrable function. "Func- 
tional calculus" can be used to show the existence and uniqueness of solutions 
of the time dependent Schrodinger and wave equation, given appropriate ini- 
tial data. Prom this point of view the time independent Green function is the 
primary object, which is unique because there is a unique selfadjoint operator 
— d^/dx^ -|- V{x) on the real line. 
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